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1  Foreward 


This  report  covers  the  period  August  15,  1995  to  December  31,  1996.  This  was  the  duration 
of  grant  DAAH04-95- 1-0527.  During  this  time,  extensive  theoretical  studies  of  hot  electron  and 
quantum  magnetotransport  in  quantum  wires  was  undertaken  under  the  author’s  supervision.  One 
graduate  student  was  supported  for  part  his  program  at  the  University  of  Notre  Dame.  He  received 
a  Ph.D.  degree  in  Electrical  Engineering  in  December  1996  and  is  currently  a  post-doctoral  research 
associate  at  the  University  of  Nebraska.  He  will  be  joining  the  Department  of  Electrical  Engineering 
at  the  University  of  California-Los  Angeles  as  a  research  associate  starting  August  1,  1997. 

This  research  resulted  in  five  journal/book  publications  and  six  conference  publications  which 
are  attached  in  the  appendix. 

A  Table  of  Contents  is  not  included  since  the  main  body  of  the  report  is  less  than  10  pages. 


2  Technical  Report 

2.1  Statement  of  problem 

Two  distint  problems  were  studied:  (a)  hot  electron  magneto-transport,  and  (b)  quantum  magneto¬ 
optical  properties  of  quantum  wires.  In  the  first  case,  we  started  with  the  Schrodinger  equation 
in  a  quantum  wire  subjected  to  a  magnetic  field.  This  equation  was  solved  numerically  using  a 
finite  difference  scheme  developed  by  us  [1]  to  yield  the  electron  wave  functions,  energy-dispersion 
relations  and  the  density  of  states.  From  these,  electron-phonon  and  electron-impurity  scattering 
rates  were  calculated  using  Fermi’s  Golden  Rule  [2-8].  These  were  then  used  in  a  Monte  Carlo 
simulator  to  find  electron  velocity  versus  field  characteristics,  electron  energy  and  other  transport 
parameters  [6,7]. 

To  study  magneto-optical  properties,  we  solved  the  Schrodinger  equation  for  an  exciton  in  a 
quantum  wire  subjected  to  a  magnetic  field.  This  was  based  on  a  variational  approach  [9-28].  The 
solutions  yield  exciton  binding  energies  and  other  parameters,  which  are  then  used  to  calculate 
bi-exciton  parameters  [10,12,14].  Third-order  dielectric  susceptibility  was  calculated  using  the 
rotating  wave  approximation.  This  third-order  non-linearity  arises  from  exciton-exciton  interaction 
and  population  saturation  of  the  excitonic  state  [11,16-18].  Furthermore,  second-order  susceptibility 
X^^^  was  also  calculated.  It  arises  from  inter-subband  transitions  and  giant  dipoles  associated  with 
them  [13,19,22,23,27].  Finally,  exciton-polariton  transport  in  an  array  of  quantum  wires  has  also 
been  investigated  [20,28]. 

2.2  Summary  of  most  important  results 
The  following  findings  were  the  most  remarkable: 

1.  It  was  found  that  a  magnetic  field  quenches  intra-subband  electron-bulk  acoustic  phonon 
scattering  rates  in  a  quantum  wire  by  several  orders  of  magnitude.  This  has  important 
implications  for  the  celebrated  integral  Quantum  Hall  effect  and  may  partially  explain  its 
robustness.  The  quenching  occurs  because  the  dominant  scattering  mechanism  is  backscat- 
tering  which  involves  transitions  between  oppositely  travelling  states.  In  the  presence  of  a 
magnetic  field,  these  states  are  localized  along  opposite  edges  of  the  wire  ( “edge  states” )  with 
very  little  overlap  between  their  wave  functions.  As  a  result,  the  matrix  element  for  transition 
between  these  states  is  drastically  reduced. 

2.  A  magnetic  field  increases  electron-surface  optical  phonon  scattering  since  the  field  skews 
the  wave  functions  of  electrons  towards  the  surfaces  (or  edges)  of  the  quantum  wire  thus 
increasing  the  overlap  between  the  surface  phonon  modes  and  the  electron  wave  function. 

3.  A  magnetic  field  has  two  different  effects  on  electron-confined  polar  optical  phonon  scattering. 
On  the  one  hand,  it  reduces  the  overlap  between  an  electron  wave  function  and  a  phonon 
mode  with  the  same  index  because  of  the  skewing  of  the  former,  but  on  the  other  hand,  it 
opens  new  channels  of  scattering  by  breaking  the  orthogonality  between  the  electron  wave 
function  and  a  phonon  mode  with  different  index.  Usually,  the  second  effect  wins,  resulting 
in  a  increase  in  the  electron-confined  optical  phonon  scattering  rates. 

4.  The  ensemble  average  transport  lifetime  in  a  quantum  wire  can  be  negative  if  the  Fermi 
velocity  (determined  by  carrier  concentration)  is  less  than  the  sound  velocity.  Such  a  negative 
lifetime  can  cause  very  low  field  velocity  overshoot. 

5.  There  can  be  an  anomalous  cooling  effect  whereby  the  electron  temperature  can  fall  below 
the  lattice  temperature  when  an  electric  field  is  applied  over  a  quantum  wire.  The  electrons 


gain  energy  from  the  electric  field,  but  lose  energy  at  a  faster  rate  by  emitting  phonons.  As 
a  result,  the  electron  gas  cools  in  an  electric  field. 

6.  The  electron  temperature  at  a  lattice  temperature  of  300  K  can  be  less  than  the  electron 
temperature  at  a  lattice  temperature  of  77  K  when  a  magnetic  field  is  present  since  phonon 
emission  (and  hence  energy  relaxation)  is  more  efficient  at  300  K. 

7.  Both  exciton  and  biexciton  binding  energies  in  a  quantum  wire  are  enhanced  by  a  magnetic 
field.  This  effect  is  opposite  to  that  caused  by  an  electric  field.  Whereas  an  electric  field 
reduces  the  binding  energy  by  exerting  equal  and  opposite  forces  on  an  electron  and  hole,  a 
magnetic  increases  the  binding  energy  by  squeezing  an  electron  and  hole  wave  function  closer 
together.  There  is  a  whole  bevy  of  magneto-optic  effects  that  are  analogs  of  well-known 
electro-optic  effects. 

8.  Both  second-  and  third-order  susceptibiltiesof  a  quantum  wire  are  enhanced  by  a  magnetic 
field.  The  magnetic  field  allows  one  to  tune  the  differential  non-linear  refractive  index  and 
absorption  in  a  quantum  wire  thereby  allowing  flexible  device  design. 

9.  A  magnetic  field  can  be  used  to  spectrally  separate  regions  of  high  refractive  index  and  high 
absorption.  As  a  result,  it  is  possible  to  sustain  polariton  transport  in  an  array  of  quantum 
wires  where  both  efficient  waveguiding  and  low-loss  propagation  is  possible.  This  has  serious 
implications  for  optical  communications. 

.3  Important  publications  resulting  from  the  grant 

1.  N.  Telang  and  S.  Bandyopadhyay,  “Magnetic  Field  Dependence  of  Energy  and  Momentum 
Relaxation  Rates  of  Hot  Carriers  in  a  Quantum  Wire  Subjected  to  a  Magnetic  Field”,  Phys. 
Low  Dimensional  Struct.  9/10,  63-74  (1996). 

2.  A.  Svizhenko,  A.  Balandin  and  S.  Bandyopadhyay  “Giant  Dipole  Effect  and  Second  Harmonic 
Generation  in  Quantum  Wires  Biased  with  a  Magnetic  Field”,  J.  Appl.  Phys.  (in  press). 

3.  A.  Balandin  and  S.  Bandyopadhyay,  “Dielectric  Confinement  of  Magneto-  Excitons  and  Biex- 
citons  in  Quantum  Wires”,  submitted  to  Zh.  Eksp.  Teor.  Fiz.:  JETP. 

4.  A.  Balandin,  A.  Svizhenko  and  S.  Bandyopadhyay,  “Non-linear  optical  mixing  in  quantum 
wires”,  to  appear  in  Quantum  Confinement  IV:  Proc.  of  the  191st  Meeting  of  the  Elec¬ 
trochemical  Society;  Fourth  International  Symposium  on  Quantum  Confinement:  Nanoscale 
Materials,  Devices  and  Systems,  Montreal,  Canada,  May  4-9,  1997. 

5.  A.  Balandin  and  S.  Bandyopadhyay,  “Polariton  transport  in  semiconductor  quantum  wires”, 
to  appear  in  Quantum  Confinement  IV:  Proc.  of  the  191st  Meeting  of  the  Electrochemical 
Society;  Fourth  International  Symposium  on  Quantum  Confinement:  Nanoscale  Materials, 
Devices  and  Systems,  Montreal,  Canada,  May  4-9,  1997. 

6.  A.  Balandin,  A.  Svizhenko  and  S.  Bandyopadhyay,  “Refractive  Index  of  a  Quantum  Wire 
Around  a  Polariton  Resonance”,  to  be  presented  at  the  Tenth  International  Conference  on 
Superlattices,  Microstructures  and  Microdevices,  Lincoln,  Nebraska,  USA. 

7.  A.  Svizhenko,  A.  Balandin  and  S.  Bandyopadhyay,  “Linear  and  Non-linear  Magneto-optical 
Properties  of  a  Quantum  Wire  Associated  with  Inter-Magnetoelectric  Subband  Resonances”, 
to  be  presented  at  the  Tenth  International  Conference  on  Superlattices,  Microstructures  and 
Microdevices,  Lincoln,  Nebraska,  USA. 


8.  A.  Balandin,  A.  Svizhenko  and  S.  Bandyopadhyay,  “Effective  Breaking  of  Inversion  Symmetry 
in  a  Quantum  Wire  with  a  Magnetic  Field”,  to  be  presented  at  the  International  Conference 
on  Nanostructures,  NANO’97,  Repino,  Russia,  June  (1997). 

9.  S.  Bandyopadhyay,  A.  Balandin  and  A.  Svizhenko,  “Dielectric  Enahancement  of  the  Oscillator 
Strengths  and  Binding  Engeries  of  Magneto-  Excitons  and  Biexcitons  in  a  Quantum  Wire” 
Spring  Meeting  of  the  American  Physical  Society,  Kansas  City,  MO,  1997. 

10.  A.  Svizhenko,  A.  Balandin  and  S.  Bandyopadhyay,  “Giant  Dipole  Effect  and  Second  Harmonic 
Generation  in  Magnetic  Field  Biased  Quantum  Wires”,  Spring  Meeting  of  the  American 
Physical  Society,  Kansas  City,  MO,  1997. 

11.  S.  Balandin  and  S.  Bandyopadhyay,  “Exciton-Polariton  Transport  in  Quantum  Wires”,  March 
Meeting  of  the  American  Physical  Society,  Kansas  City,  MO,  1997. 

2.4  Personnel 

1.  Prof.  Supriyo  Bandyopadhyay,  principal  investigator 

2.  Dr.  Nina  Telang,  Research  Engineer,  Motorola,  Inc,  Austin,  TX 

3.  Dr.  Alexander  Balandin,  post-doctoral  associate.  University  of  Nebraska 

Dr.  Balandin  was  a  graduate  student  working  on  this  grant  who  received  his  Ph.D.  in  December 
1996. 
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the  giant  dipole  effect  in  magnetic-Iield-biased  semiconductor  quantum 
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electric  subbands  within  the  conduction  band;  some  of  these  transitions 
are  forbidden  in  the  absence  of  the  magnetic  field.  The  possibility  of  sec¬ 
ond  harmonic  generation  in  a  quantum  wire  biased  with  a  magnetic  field 
has  also  been  examined.  We  will  show  that  the  simultaneous  presence  of 
a  symmetric  confinement  (electrostatic)  potential  and  an  external  mag¬ 
netic  field  may  lead  to  uneven  charge  distribution  along  the  width  of  the 
wire  which  breaks  the  inversion  symmetry  in  a  generic  quantum  well  or 
wire.  This  results  in  a  non-zero  third-order  susceptibility  and  associated 
optical  nonlinearity.  This  work  was  supported  by  the  US  Army  Research 
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ALEXANDER  BALANDIN,  Department  of  Electrical  Engineering,  Uni¬ 
versity  of  Nebraska,  Lincoln,  NE.,  ALEXEI  SVIZHENKO,  Department 
of  Electrical  Engineering,  University  of  Notre  Dame,  Notre  Dame,  IN. 

—  We  will  report  a  partly  variational  calculation  of  magneto-  exci- 
ton  and  biexciton  (excitonic  molecule)  binding  energies  in  a  rectangular 
semiconductor  quantum  wire  taking  into  account  dielectric  confinement 
effects.  The  discontinuity  of  the  dielectric  constant  at  the  boundaries  be¬ 
tween  the  wire  and  surrounding  medium  causes  redistribution  of  electric 
field  which  leads  to  higher  binding  energies  of  excitons  and  biexcitons. 

The  dielectric  confinement  effects  are  introduced  via  the  image-charge 
method.  Strong  dielectric  enhancement  of  the  exciton  binding  energies 
(up  to  2.5  times)  and  oscillator  strengths  is  found  for  quantum  wires  of 
the  dimensions  comparable  to  the  exciton  Bohr  radius.  The  enhance¬ 
ment  decreases  substantially  in  a  magnetic  field  because  of  a  charge 
cancellation  effect.  This  work  was  supported  by  the  US  Army  Research 
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Exciton-Polariton  Transport  in  an  Array  of  Semicon¬ 
ductor  Quantum  Wires  ALEXANDER  BALANDIN,  SUPRIYO 
BANDYOPADHYAY,  Department  of  Electrical  Engineering,  University 
of  Nebraska,  Lincoln,  NE.  —  Polariton  effects  play  an  important  role 
in  low  temperature  absorption  of  bulk  and  nanostructured  semiconduc¬ 
tors  in  the  spectral  range  close  to  the  exciton  resonance  frequencies. 
The  formation  of  polaritons  (exciton-photon  coupled  states)  modifies 
the  transport  of  light  through  the  medium.  In  particular,  the  medium 
of  propagation  becomes  substantially  more  transparent  provided  that 
the  exciton-polariton  coherence  length,  critical  temperature  and  exciton 
decay  parameter  meet  certain  criteria.  We  wiU  show  that  in  a  properly 
designed  array  of  quasi  one- dimensional  structures  (quantum  wires)  of 
dimensions  comparable  to  the  exciton  Bohr  radius,  the  application  of 
an  external  magnetic  field  can  strongly  increase  the  critical  exciton- 
polariton  coherence  length  and  the  critical  temperature.  Our  numerical 
results  are  based  on  a  partially  variational  calculation  of  the  exciton 
longitudinal-transverse  splitting  as  a  function  of  wire  dimensions  and 
magnetic  flux  density.  The  temperature  dependence  is  evaluated  un¬ 
der  the  assumption  that  interactions  with  acoustic  and  polar  optical 
phonons  are  the  dominant  scattering  mechanisms  for  excitons.  This 
work  was  supported  by  the  US  Army  Research  Office  under  contract 
DAAH04-95- 1-0586  and  DAAH04-95- 1-0527. 
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Dielectric  Confinement  of  Magneto-  Excitons 
and  Biexcitons  in  Quantum  Wires 

A.  Balandin*and  S.  Bandyopadhyay 
Department  of  Electrical  Engineering 
University  of  Nebraska 
Lincoln,  Nebraska  68588  U.S.A. 


The’-  effect  of  dielectric  confinement  on  magneto-  excitons  and  biexcitons 
in  a  quantum  -wire,  surrounded  by  a  medium  -with  lo-wer  dielectric  cosntant, 
is  investigated  theoretically.  For  free-standing  wires,  very  strong  dielectric 
enhancement  of  binding  energies  and  oscillator  strengths,  along  with  a  con¬ 
comitant  decrease  in  radii,  is  found  for  wires  whose  widths  are  comparable  to 
the  effective  Bohr  radius  of  the  exciton.  In  GaAs  wires,  the  combination  of 
dielectric  confinement,  quantum  confinement  and  magnetostatic  confinement 
can  make  both  the  exciton  and  biexciton  binding  energies  exceed  the  room- 
temperature  thermal  energy  for  realistic  wire  widths.  We  also  show  how  a 
magnetic  field  can  be  used  to  experimentally  separate  the  effect  of  dielectric 

confinement  and  quantum  confinement. 
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I.  INTRODUCTION 


Quasi  one-dimensional  structures  (quantum  wires)  exhibit  enhanced  ex- 
citon  and  biexciton  binding  energies  and  oscillator  strengths  which  endow 
them  with  strong  non-linear  optical  properties^"®.  The  enhancements  ac¬ 
crue  primarily  from  quantum  confinement  which  is  spatial  confinement  of 
the  electron  and  hole  wavefunctions  in  a  quasi  one-dimensional  region.  This 
confinement  can  be  increased  further  by  a  magnetic  field  which  “squeezes” 
the  electron  and  hole  wave  functions  into  edge-states  or  cyclotron  orbits.  At 
the  same  time,  if  the  medium  surrounding  the  wire  has  a  smaller  dielec¬ 
tric  constant,  then  the  effective  dielectric  constant  of  the  entire  system  is 
reduced  which  then  reduces  the  screening  of  the  attractive  interaction  be¬ 
tween  an  electron  and  a  hole.  As  a  result,  an  electron  and  a  hole  become 
even  more  tightly  bound  and  this  increases  the  binding  energy  and  oscilla¬ 
tor  strengths  further®"®.  In  this  paper,  we  examine  this  latter  effect  -  the 
so-called  “dielectric  confinement”  -  in  the  presence  of  a  magnetic  field. 

While  dielectric  confinement  effects  in  quantum  wells  and  dots  have  been 
studied®’^®,  scant  attention  has  been  given  to  quantum  wires.  So  far,  the  very 
few  theoretical  treatements  that  have  been  reported  for  wires  have  concen¬ 
trated  narrowly  on  specific  geometries  and  carrier  density  regimes.  In  ref. 
[11],  the  effect  of  dielectric  confinement  on  exciton  binding  energy  was  found 
variationally  for  cylindrical  wires  whose  radii  R  were  much  smaller  that  the 
exciton  Bohr  radius  as-  The  mathematical  model  in  this  work  is  valid  only 
when  R  «  as  and  is  therefore  of  limited  use.  Ref.  [12]  treated  a  wedge 
shaped  wire.  We  consider  a  wire  of  rectangular  cross-section.  To  our  knowl¬ 
edge,  ours  is  also  the  first  study  where  the  effect  of  an  external  magnetic 
field  has  been  taken  into  account.  Furthermore,  we  have  also  examined  the 
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effect  of  dielectric  confinement  on  biexcitons.  This  is  an  extremely  important 
topic  since  exciton-exciton  interaction  leading  to  the  formation  of  biexcitons 
(and  possibly  higher  order  excitonic  complexes)  is  the  major  source  of  opti¬ 
cal  non-linearity  in  quantum  wires^^’^^.  Increased  biexciton  binding  energies 
translate  into  increased  third-order  dielectric  susceptibility  (and  therefore 
increased  optical  non-linearity)  in  quantum  wires.  Finally,  any  increase  in 
biexciton  binding  energy  makes  these  complexes  more  stable  and  therefore 
improves  the  chances  of  experimental  observation. 

This  paper  is  organized  as  follows.  In  the  next  section,  we  present  the 
theory.  This  is  followed  by  results  and  finally  the  conclusion. 

II.  THEORY 

A.  Dielectrically  confined  magneto-excitons 

We  consider  a  quantum  wire  of  rectangular  cross-section  as  shown  in  Fig. 
1  (a)  with  infinite  potential  barriers  located  at  y  =  '^Ly/2  and  z  —  ±1/^/2. 
A  magnetic  field  of  flux  density  B  is  applied  along  the  ^-direction.  For 
nondegenerate  and  isotropic  bands,  the  Hamiltonian  of  a  free  Wannier-Mott 
exciton  in  this  system  is  given  within  the  envelope-function  approximation 
by 

_  Px  .  Vx  ,  VyH^fz^ 

+  +  eB{ye/me  -f-  yhfmh)px  +  {vl/me  +  yl/ rrih) 

"b  Uc(^Xe,  Xfi,  yei  Vtii  ^h)  "b  Vh^  ^h)  "b  yh,  ^h)y 


where  we  have  chosen  the  Landau  gauge  for  the  magnetic  vector  potential 


The  quantities  m^,mh,  {xe,h,ye,h,Ze,h)  are  the  effective  masses  (coordinates) 
of  electrons  and  holes  respectively,  1/Ai(=  l/nie  +  l/mh,)  is  the  reduced  mass, 
M{=me  +  mh)  is  the  total  mass,  Uc{xe,  Xh,  Ve,  Vh,  Ze,  is  the  electron-hole 
Coulomb  interaction  energy,  Us{ye,yh,Ze,Zh)  is  the  spatial  confinement  po¬ 
tential  energies  for  electrons  and  holes  along  the  transverse  y  and  z  direc¬ 
tions,  and  UD{ye,yh,Ze,Zh)  is  the  dielectric  confinement  for  electrons  and 
holes  along  y  and  z  directions  arising  from  the  redistribution  of  the  electric 
field  due  to  the  discontinuity  of  the  dielectric  constant  at  the  wire  inter¬ 
faces.  We  neglect  Zeeman  splitting  since  the  Lande  g-factor  is  small  in  most 
technologically  important  semiconductors. 

For  convenience,  we  replaced  Xg./i-coordinates  in  Eq.  (1)  by  the  center- 
of-mass  {X)  and  relative  coordinate  (x): 

X  =  {meXe+mhXh.)/M, 

X  =  Xe-  Xh. 

For  free-standing  quantum  wires  (surrounded  by  air  or  vacuum),  we  can 
assume  hard-wall  boundary  conditions: 

Usiye^yhi  Zg,  Z}i)  =  0,  Ze,h  ^  -^z  A  ye,h  ^  Ly  (2) 

=  oo,  otherwise. 

Finite  barrier  heights  introduce  only  small  corrections  for  wires  with  width 
and  thickness  >  50  A;  therefore,  our  results  are  quite  general  for  wires  of 
realistic  widths. 

The  Coulomb  interaction  between  any  two  charged  particles  in  a  quantum 
wire  is  given  by  its  full  three-dimensional  form 

Uc{Xg,XH,ye,yh,Z,,Zh)  =  _  2.^)2  +  > 

(3) 
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where  Ciu  is  the  dielectric  constant  of  the  wire. 

To  calculate  the  dielectric  confinement  energy,  Ud,  we  adopt  the  image- 
charge  method  which  is  known  to  provide  good  qualitative  agreement  with 
experiments®’^®’^®.  The  image-charge  method,  which  is  well-established  in 
electrostatics,  accounts  for  the  electric  field  induced  by  charged  particles  in 
parallel  or  spherical  geometries  by  means  of  imaginary  charges  placed  in  sur¬ 
rounding  media.  The  values  of  the  image-charges  are  determined  from  the 
continuity  conditions  for  the  electrostatic  potentials  and  the  normal  compo¬ 
nents  of  the  displacement  vector  at  the  interfaces.  Following  the  derivation 
of  Refs.  [9,  10],  we  can  write  the  dielectric  confinement  potential  of  an 
electron-hole  pair  within  the  quantum  wire  as 


Uo  =  ul%  +  U%  +  ui7,'‘\ 


(4) 


where 


p2  00  00 


“  2e 


U 


J=-oom=-oo  [{Ve  ye,l,w)^  + 

2e,j,  i-^QQm=-<x  [(2/^  ~  yh,l,rnf'  +  i^h  ~~ 


(5) 

(6) 


U, 


(e-h)  _  _£ 
attr 


,2  00  00 


=  -fE  E 


i  “oomfl^oo  [i.^e  “  +  (^/e  “  yh,l,mY  +  (^e  “  Zh,l,mYYI'^  ' 


(7) 


and 


2^0!,/,m)  —  (^^aj  Tj/ -k  (  1)  L^TH  (  1)  ZfY), 


with  a  =  eov  h]  and  l,m  —  — oo, ...,  oo.  The  l  =  m  =  0  term  is  excluded  from 
the  summation  in  Eqs.  (5)-(7).  The  dielectric  constant  of  the  surrounding 
medium  is  assumed  to  be  less  than  the  dielectric  constant  of  the  wire  material 
(cfr  <  Cty),  which  is  usually  the  case  in  real  systems.  Finally,  the  quantity 
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^  {^w  —  ^b)/{'^w  +  £&)  is  a  measure  of  the  dielectric  misfit  between  the 

quantum  wire  material  (dielectric  constant  e^)  and  the  boundary  material 
(dielectric  constant  ef). 

One  can  see  from  Eq.  (4)  that  the  dielectric  confinement  term,  C/d,  is 
written  as  a  sum  of  two  positive  self-energy  terms  resulting  from  repulsive 
interactions  of  electron  -“electron  image”  (hole  -  “hole  image”);  and  a  neg¬ 
ative  term  which  is  due  to  attractive  interaction  between  electron  -  “hole 
image”  (hole  -  “electron  image”).  The  summations  over  the  infinite  series  in 
Eqs.  (5)  -  (7)  come  about  because  any  charge  image,  in  turn,  creates  another 
image  in  the  opposite  flat  boundary  of  the  wire  (see  Fig.  1  (b)). 

For  materials  with  a  small  dielectric  misfit  (e.g.  CdS  dispersed  in  alumina 
for  which  ^  0.5),  one  can  retain  only  the  first  order  terms  {\l\ -f  |ml  =  1)  in 

Eqs.  (5)  -  (7);  for  material  systems  with  higher  dielectric  misfit  (such  as  GaAs 
surrounded  by  air  for  which  ^  0.9),  higher  order  terms  are  required.  For  the 

sake  of  illustration,  we  will  write  the  first  order  terms  in  C/d  explicitly.  Setting 
I  =  0  and  m  =  ±1  we  obtain  the  terms  which  correspond  to  interactions  of 
the  electron  and  hole  with  their  images  in  the  flat  boundaries  along  z  direction 
(Lz  is  the  thickness  of  the  wire).  Similarly,  setting  m  =  0  and  /  =  ±1  we 
obtain  the  terms  which  correspond  to  interactions  of  the  electron  and  hole 
with  their  images  in  the  flat  boundaries  along  y  direction  {Ly  is  the  width  of 
the  wire) .  Combining  these  terms  together,  we  get 

7-/(e)  _  1  —  gfew  1  1 

^2e^\^  +  eb’^\2z,-Lz\  \2z,  +  Lz\ 

1  1 
\2^ye- Ly\'^  \2ye  +  Ly\’ 


TjW 

^self 


—  _ f*)/ _ - _ 1 _ I - 

2€yj  C’uj  Cb  \2zii  -k  Lz  \2zh  —  Lz 


(9) 
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The  exciton  ground  state  in  the  system  described  by  the  Hamiltonian 
of  Eq.  (1)  is  found  variationally  by  minimizing  the  expectation  value 
<  >  where  the  exciton  wave  function  is  properly  normalized. 

We  will  be  interested  only  in  excitonic  states  that  can  be  accessed  optically. 
Under  optical  excitation,  the  center-of-mass  motion  can  be  neglected  since 
the  photon  momentum  is  too  small  to  create  states  of  significant  center-of- 
mass  kinetic  energy.  Therefore,  the  wave  function  of  an  optically  generated 
exciton  (s  symmetry  only)  with  center-of-mass  momentum  =  0  is  given 
by 

^  Ve,  Vh,  Ze,  Zh)  =  gt{z:,  T])'tpe{ye,  Ze)'lphiyh,  Zh)  (11) 

=  9t{x,  V)(f>eiye)Myh)Xe{Ze)Xh{Zh), 

where  g  is  the  x-component  of  the  exciton  wave  function,  and  (p  and  x  are 
respectively  the  y-  and  z-components  of  independently  confined  electron  and 
hole  making  up  an  exciton. 

The  separation  of  the  wave  function  into  a  product  of  g,  (p  and  x 
Equation  (11)  is  permissible  only  in  the  limit  of  weak  Coulomb  interaction 
between  the  charged  particles,  i.e.  for  materials  with  relatively  large  dielectric 
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constants  (narrow  gap  semiconductors).  The  z-component  Xe,h{ze,h)  is  not 
affected  by  a  magnetic  field  oriented  along  the  z-direction.  Therefore,  it  is  a 
particle-in-a-box  state: 

Xe,h{z^,h)  =  (12) 

The  y-component  of  the  electron  and  hole  wave  functions,  (j)e,h{ye,h),  are  af¬ 
fected  by  the  magnetic  field  and  are  calculated  numerically.  This  is  done 
by  solving  the  one-particle  Schrodinger  equation  using  a  finite  difference 
scheme.^® 

Finally,  the  x-component  of  the  wave  function  (i.e.  the  component  along 
the  free  direction  or  wire  axis)  is  given  by  a  Gaussian-type  “orbital”  function: 

=  (13) 

The  ground  state  exciton  binding  energies  can  now  be  found  using 
the  relation® 


+  BSi,  +  E%  +  (14) 

Here  E^j  are  the  lowest  electron  and  the  highest  heavy  hole  magneto¬ 
electric  subband  energies  in  the  quantum  wire  measured  from  the  bottom  of 
the  bulk  conduction  band  and  the  top  of  the  bulk  valence  band  respectively. 
The  total  energy  E*°^  is  determined  by  minimizing  the  expectation  value  of 
the  total  Hamiltonian  by  varying  the  parameter  t],  i.e., 

E^°*  =  mirir,  <  >,  (15) 

and  the  electron  self-energy  E^fif  is  given  as 

^self  ^e(ye)  ■2(e)|^ieZ/lV’e(l/e>  ■2'e)  ^  •  (16) 
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The  hole  self-energy  is  obtained  from  the  last  expression  by  substituting  the 
hole  index  “h”  in  place  of  the  electron  index  “e.” 

The  evaluation  of  self-energy  integrals  in  Equation  (16)  is  not  straight¬ 
forward  because  of  the  unphysical  divergences  associated  with  the  “poles” 
of  the  self-energies  in  Equations  (8)  and  (9).  These  divergences  are  a  conse¬ 
quence  of  assuming  that  the  dielectric  constant  changes  abruptly  at  the  wire 
interface  (this  is  also  equivalent  to  assuming  that  the  induced  surface  charge 
has  no  spatial  extension).  In  order  to  avoid  these  divergences,  we  adopt  the 
approach  of  ref.  [19]  and  assume  that  the  surface  charge  is  located  symmet¬ 
rically  with  respect  to  the  boundary  and  has  a  spatial  extension  5  less  than 
the  lattice  constant.  The  self-energy  terms  are  then  calculated  as  “principal 
values”  with  the  replacement 

__i _ .  \2z.±LA 

\2z,±L,\  (2z,  ±  L,)2  +  52 

Unlike  the  self-energy  integrals,  the  Coulomb  integrals  <  >  and 

<  ^lUattrl^  >  are  “proper”  integrals  as  long  as  Uc  and  Uaur  of  Eqs.  (3)  and 

(10)  are  treated  in  their  full  three-dimensional  form.  Despite  the  character¬ 
istic  1/r  singularity,  these  integrals  do  not  diverge  because  in  two  or  three 
dimensions,  the  surface  element  2Trrdr  or  the  volume  element  goes 

to  zero  as  rapidly  as,  or  more  rapidly  than,  r.  Consequently,  the  integrals 

<  >  and  <  $|17aHr|^  >  are  always  finite  if  Uc  and  Uattr  are  treated 
in  their  three-dimensional  form.  We  cannot  do  the  same  for  the  self-energy 
terms  because  they  depend  only  on  one  coordinate  -  z  ov  y  -  and  hence  the 
integration  is  effectively  one-dimensional.  Therefore,  it  becomes  necessary 
to  introduce  the  artificial  parameter  5.  Nonetheless,  we  always  verify  that 
the  final  result  of  any  calculation  does  not  depend  sensitively  on  S.  Typi¬ 
cally,  several  orders  of  magnitude  variation  in  6  introduces  only  a  few  percent 
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variation  in  the  final  result. 


B.  Dielectrically  confined  magneto-biexcitons 


In  order  to  study  the  effect  of  dielectric  confinement  on  magneto- 
biexcitons,  we  use  the  same  quantum  wire  geometry  as  in  Fig.  1(a).  Our 
mathematical  model  will  be  based  on  several  assumptions.  First,  we  consider 


the  regime  of  intermediate  carrier  densities,  n,  for  which  the  Wigner-Seitz 
radius  ps  =  (3/47rn)^/®  »  cb,  where  is  again  the  Bohr  radius  of  the 
exciton.  Extremely  high  carrier  densities  (excitation  levels)  may  lead  to  dis¬ 
sociation  of  all  excitonic  complexes  and  formation  of  the  electron-hole  plasma 


which  is  not  within  the  scope  of  this  work.  Second,  we  neglect  any  polariton 
effects  which,  in  fact,  may  be  important  for  very  narrow  wires.^°’^^ 

For  nondegenerate  and  isotropic  bands,  the  Hamiltonian  of  a  biexciton  in 


a  quantum  wire  subjected  to  a  magnetic  field  is  given  in  the  envelope-function 
approximation  by^^ 

52 


H^x  = 


2/i  dx\^ 


-b 


dxlb 

52 


n-  0“  0“  a-  . 

2me  dyi^ 


5^  . 

nih^dxlf^  dXiadXab  dX2bdXab 


+ 


92  02 

+ 


92 

+  — 7;) 


2mhdya^  dyt^  dza^  '  dz^^' 
e^B^^yl  +  yj  yl  +  Vb^  eBih  92 
2  m.1.  m.  ^^x^. 


92 

2  ^  me  '  mh  '  '  TUe  ''"'■9a;ia  ”*"^^9x26 

eBin,  92  ,  -  92  ,  92  , 


rUh  ''''  dXab 

-b  Vc(Xe,  Xh,  ye,  yh,  Ze,  Zh)  +  V's(y, 


“9X26' 

rt  ^  4-  (m 


(18) 


where  Vs  is  the  spatial  confinement  potential  along  the  y  and  z  directions, 
Vc  is  the  Coulomb  interaction  between  various  charged  entities,  and  Vd  is 
the  dielectric  confinement  for  all  four  particles  comprising  a  biexciton.  The 
electron  coordinates  bear  numerical  subscripts  .1  and  2  while  the  hole  co¬ 
ordinates  bear  alphabetical  subscripts  a  and  b.  For  convenience,  we  again 


11 


replace  the  coordinates  in  Eq.  (18)  by  the  center-of-mass  (X)  and  relative 
coordinates  {Xa,-y)  using  the  following  transformations'^’^® 

X  =  [Tne{xi+X2) +  mh{xa  +  Xb)]/2{me  +  mh),  (19) 

^ICL  ^1  ^CL 

^26  =  3:2—  Xb 

XqI)  Xq  Xfjj 


where  a,  7  =  1, 2,  a,  b,  and  a  ^  7. 

As  in  Eq.  (2),  the  spatial  confinement  potential  Vs  is  given  by 

Vsiyi,a,y2,b,Zi,a,Z2,b)  =  0,  Zi^a,2,b  <  yi,a,2,b  <  Ly  (20) 

=  oo,  otherwise. 


The  Coulomb  interaction  term  is  given  by 


^Coulomb  —  i  ■{ 


.  (-1)  .  (-1)  ,  (-1) 


+ 


+ 


yj^ia  +  rh  y/^ib  +  ^ib  yj^lb  +  rjb 

(-1)  1  1  1 
^  ^  +-7--- . --+-7-- .  }, 


+  \lAb  +  'rlb 


where 


(21) 


I'i,  =  (Vt.  -  Vjf  +  (*0  -  hf: 


Xa,^  —  Xa  Xy, 

with  a,  7  =  1, 2,  a,  &,  and 

The  dielectric  confinement  term  Vd  is  analogous  to  Eqs.  (4) -(7).  The 
only  principally  different  terms  which  appear  in  the  biexciton  problem  and 
are  not  present  in  the  exciton  problem  are  the  repulsive  electron  -  “another 
electron-image”  (hole  -  “another  hole-image”)  interactions.  These  terms  have 
the  same  structure  as  the  attractive  terms  but  bear  the  opposite  sign.  For 
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the  sake  of  illustration,  we  have  written  the  first-order  terms  explicitly  in  the 
Appendix. 


The  variational  procedure  for  calculating  the  biexciton  binding  energies 
is  the  same  as  that  reported  in  ref.  [22,  23].  The  trial  wave  function  is  chosen 
to  be  a  singlet  state  with  the  electron-hole  pair  contributions  given  by  the 
Gaussian-type  “orbital”  functions: 

^  =  ;^{^lo^26  +  '>p2ai’lb}Gab{Xab),  (22) 

where 

'0Q7  ~  9a'y^a{ya)Xa{^a)^y{y'i/)X‘y{^'y) t  (23) 

with 


Gab{Xab)^e-^l’>/-\  (24) 

and 

gia  =  e-</^\  (25) 

=  (26) 

=  (27) 

=  (28) 

The  wave  functions  ga-y  are  Gaussian  orbitals  whose  “spread”  rj  and  r  are 
variational  parameters.  It  is  obvious  that  these  parameters  physically  cor¬ 
respond  to  the  electron-hole  and  hole-hole  separations  along  the  length  of 
the  quantum  wire.  The  quantity  S'  is  a  normalization  constant  and  ^e,h{ye,h) 
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[Xe,h{ze,h)\  y-componeiits  [z-components]  of  the  wave  functions  of  inde¬ 
pendently  confined  electrons  and  holes  as  in  the  previous  section. 

In  our  own  previous  calculation  of  the  magneto-biexciton  binding  energies 
in  ref.  [22],  we  employed  two  variational  parameters;  77  and  r.  The  use  of  only 
two  variational  parameters  tacitly  assumes  that  the  electron-hole  separation 
in  a  biexciton  is  the  same  as  that  in  a  free  exciton.  In  the  present  work, 
we  do  away  with  this  assumption  and  improve  on  our  previous  calculations 
by  including  a  third  variational  parameter  rf,  which  allows  the  electron-hole 
pair  within  the  biexciton  to  relax.  This  is  done  following  Ref.  [24].  The 
biexciton  binding  energy  is  found  from  the  relation 

>  -2mm,,  <  ip\H^\ip  >,  (29) 

where  minr,r^»  <  >  is  found  by  minimizing  the  expectation  value 

of  the  Hamiltonian  with  respect  to  the  hole-hole  variational  parameter 

r  and  the  electron-hole  relaxation  parameter  rf,  while  mirir,  <  >  is 

found  as  described  in  the  previous  section. 

Before  concluding  this  section,  we  should  make  a  few  remarks  about  the 
numerical  details.  As  we  have  already  mentioned  in  the  previous  section, 
the  Coulomb  terms  treated  in  their  full  three-dimensional  form  do  not  cause 
divergence,  and  can  be  handled  without  introducing  any  ad  hoc  truncation 
procedure.  However,  this  strategy  is  also  CPU-costly  since  the  number  of 
independent  coordinates  (spatial  dimensions)  increases  threefold.  The  cal¬ 
culation  of  the  expectation  value  of  the  biexciton  Hamiltonian  involves  eval¬ 
uation  of  multi-dimensional  integrals  (n  =  9)  in  real  space.  The  evaluation 
is  achieved  by  iterated  applications  of  product  Gauss  formulas.  The  integral 
is  first  estimated  by  a  two-point  tensor  product  formula  in  each  coordinate 
direction.  Then  for  i  =  1,  ...,n  the  routine  calculates  a  new  estimate  by  dou- 
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bling  the  number  of  points  in  the  t-th  direction,  but  halving  the  number 
immediately  afterwards  if  the  new  estimate  is  not  appreciably  different. 

III.  RESULTS  AND  DISCUSSION 

A.  Magnetoexciton  and  magnetobiexciton  binding 
energies 

All  the  results  that  we  present  in  this  section  are  relevant  to  free-standing 
GaAs  quantum  wires  surrounded  by  vacuum  or  air.  Electron  and  hole  effec¬ 
tive  masses  are  assumed  to  be  0.067mo  and  0.5mo  respectively.  The  relative 
permittivity  of  GaAs,  e^,  is  taken  as  12.9.  The  dielectric  misfit  between 
GaAs  and  vacuum  is  estimated  to  be  ^  =  0.856. 

The  binding  energies  of  dielectrically  confined  magneto-excitons  are  cal¬ 
culated  from  Eqs.  (14)  -  (16).  In  order  to  give  some  idea  about  the  relative 
scale  of  the  terms  in  Eq.  (14),  we  present  in  Fig.  2  the  confinement  energies 
(the  magneto-electric  subband  bottom  energies  and  as  functions 

of  the  wire  width  for  two  different  values  of  a  magnetic  fiux  density.  The  fun¬ 
damental  band  gap  energy  Eq  in  Fig.  2  is  assumed  to  be  zero  for  clarity. 
Although  not  shown  in  the  figure,  the  values  of  the  self-energy  terms  are 
about  20-30%  higher  than  the  corresponding  confinement  energies. 

In  Fig.  3,  we  present  the  ground  state  binding  energy  of  a  dielectrically 
confined  exciton  as  a  function  of  wire  width  Ly.  For  comparison,  the  binding 
energy  without  dielectric  confinement  effects  is  shown  by  the  dashed  line. 
The  dash-dotted  curve  in  Fig.  3  represents  the  case  when  only  first  order 
terms  in  Eq.  (5)-(7)  were  taken  into  account.  As  one  can  see  the  discrepancy 
between  this  curve  and  the  solid  curve  which  is  obtained  using  higher  order 
terms  (up  to  4th)  is  quite  pronounced.  This  difference  is  always  appreciable 
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for  material  systems  with  large  values  of  dielectric  misfit  ^  >  0.5.  When 
the  dielectric  misfit  is  not  large  (e.g.  the  GaAs/AlGaAs  system  rather  than 
the  GaAs/vacuum  system),  just  the  first  order  correction  may  be  enough. 
The  binding  energy  always  decreases  with  increasing  wire  width  for  all  cases 
as  expected  because  of  decreasing  spatial  (electrostatic)  confinement.  For  a 
200  A  thick  wire,  the  exciton  ground  state  binding  energy  is  enhanced  by 
approximately  2.7  times  because  of  dielectric  confinement.  Note  that  the 
binding  energy  is  larger  than  the  room-temperature  thermal  energy  (kT  = 
26  meV)  for  wire  widths  less  than  300  A. 

The  exciton  radius  (or  more  correctly  the  longitudinal  length)  along  the 
X  direction  is  given  analytically  as  p  =  V<  >  =  77/2^/^  for  the  chosen 
trial  wave  function  in  Equation  (11).  This  radius  is  shown  in  Fig.  4  as  a 
function  of  the  wire  width  with  and  without  dielectric  confinement  taken  into 
account.  It  is  interesting  to  note  that  the  ratio  (/?  =  p  —  Pd)/p  shows  very 
weak  dependence  on  the  wire  width  Ly.  Here  po  is  the  the  exciton  radius 
in  a  system  with  dielectric  confinement  and  p  is  the  same  without  dielectric 
confinement. 

In  Fig.  5  we  present  the  dependence  of  the  exciton  radius  on  magnetic 
flux  density  for  two  cases:  (i)  with  dielectric  confinement,  and  (ii)  without 
dielectric  confinement.  The  important  thing  to  note  is  that  the  difference 
between  the  two  curves  decreases  with  increasing  magnetic  flux  density.  The 
exciton  radii  ratio  P  depends  quite  strongly  on  magnetic  field.  For  the  wire 
with  dimensions  Ly  x  =  500A  x  200A,  the  ratio  can  be  approximated  as 
P  =  0.128  —  0.0062jB,  where  B  is  in  tesla.  This  monotonically  decreasing 
dependence  of  P  on  B  can  be  easily  explained.  A  magnetic  field  squeezes 
the  exciton  wave  function  along  all  directions,^’^^  pushing  the  electron  and 
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the  hole  closer  together.  As  a  result,  the  exciton  increasingly  approximates 
a  truly  electrically  neutral  uncharged  single  particle  rather  than  a  bound 
electron-hole  pair.  Since  a  chargeless  particle  has  no  image  charge,  the  di¬ 
electric  confinement  effect  falls  off  rapidly. 

The  strong  dependence  of  /?  on  magnetic  field  has  a  very  important  im¬ 
plication.  It  appears  that  because  of  this  feature,  a  magnetic  field  can  be 
used  to  separate  the  effects  of  spatial  confinement  and  dielectric  confinement 
on  the  exciton  radius.  This  could  not  be  done  by  varying  the  wire  width  or 
spatial  confinement  (e.g.  in  split-gate  quantum  wires)  since  evidently  ^  does 
not  depend  sufficiently  strongly  on  Ly  (recall  Fig.  4). 

Fig.  6  shows  the  corresponding  dependence  of  the  exciton  binding  energy 
on  the  magnetic  flux  density  for  the  same  wire  dimensions  as  in  Fig.  5.  As  one 
can  see,  the  difference  between  binding  energies  with  and  without  dielectric 
confinement,  ,  does  not  depend  on  the  magnetic  field  as  strongly  as  the 
difference  between  the  radii.  The  magnitude  of  the  dielectric  enhancement 
of  the  binding  energy  is  in  a  good  qualitative  agreement  with  experimental 
results  published  for  quantum  wells®^^.  Unfortunately,  no  data  are  available 
yet  for  quantum  wires  so  that  a  quantitative  comparison  is  not  possible. 

The  increase  in  the  binding  energy  and  the  decrease  in  the  radius  as  a 
result  of  dielectric  confinement  can  be  simply  ascribed  to  a  lowering  of  the 
effective  dielectric  constant  (and  hence  the  screening)  in  the  entire  system 
due  to  the  image  charges.  A  magnetic  field  reduces  the  image-charge  effect 
because  it  squeezes  the  electron  and  hole  togther  thereby  tending  to  con¬ 
dense  them  into  an  uncharged  single  particle.  Cosnequently,  a  magnetic  field 
quenches  the  dielectric  confinement. 

In  Fig.  7  we  present  the  dependence  of  the  exciton  radius  on  dielectric 
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misfit  with  and  without  the  magnetic  field  present.  As  one  should  expect, 
the  radius  decreases  with  increasing  misfit  for  both  cases  provided  that  the 
dielectric  constant  of  the  surrounding  material  is  lower  than  that  of  the  wire. 
This  is  simply  a  validation  of  the  fact  that  a  higher  dielectric  misfit  results 
in  a  higher  dielectric  confinement.  Furthermore,  the  exciton  radius  is  more 
sensitive  to  dielectric  confinement  if  there  is  no  magnetic  field.  The  physics  of 
it  is  intuitively  obvious:  the  wave  function  of  the  exciton  in  a  quantum  wire 
subjected  to  the  field  (10  tesla)  is  already  compressed  due  to  strong  magneto¬ 
static  confinement  and  does  not  “feel”  the  additional  dielectric  confinement 
as  much. 

The  corresponding  dependence  of  the  exciton  binding  energy  on  the  di¬ 
electric  misfit  is  shown  in  Fig.  8.  It  is  interesting  to  note  that  the  dependence 
is  almost  exactly  linear,  for  the  magnetic  flux  density  of  10  tesla  it  can  be 
approximated  by  the  interpolation  formula  Eg  ~  9.240  d- 16.376  ^  -1-0.089 
where  the  pre-factor  of  the  quadratic  term  is  more  than  two-orders  of  mag¬ 
nitude  smaller  than  that  of  the  linear  term. 

In  Fig.  9  (a)  we  present  the  magneto-biexciton  binding  energy  as  a  func¬ 
tion  of  the  magnetic  field  for  two  cases:  with  dielectric  confinement  and 
without  it.  As  we  can  see,  the  overall  behavior  of  the  binding  energy  is 
similar  to  that  in  the  case  of  magneto-excitons,  although  the  dielectric  con¬ 
finement  efi'ect  is  not  so  pronounced.  The  weaker  effect  is  probably  caused 
by  the  fact  that  there  are  like-charged  entities  (two  holes  and  two  electrons) 
in  a  biexciton  that  repel.  As  a  result,  the  dielectric  confinement  effect  is  not 
so  strong.  Fig.  9  (b)  shows  the  biexciton  binding  energy  as  a  function  of  the 
wire  width.  It  is  interesting  to  note  that  in  the  presence  of  dielectric  confine- 
men  (“-f”  marked  curve)  the  binding  energy  exceeds  the  room-temperature 
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thermal  energy  of  26  meV  for  wire  width  less  than  100  A.  Application  of  an 
external  magnetic  field  allows  to  exceed  room-temperature  thermal  energy 
for  even  wider  wires. 


III.  DIELECTRIC  ENHANCEMENT  OF  OSCILLA¬ 
TOR  STRENGTH 


In  this  section,  we  examine  the  effect  of  dielectric  confinement  on  os¬ 
cillator  strength.  In  order  to  calculate  the  oscillator  strength  of  the  exciton 
transition,  we  need  to  evaluate  the  momentum  matrix  element  which  is  given 
by25 

=  1^/  dk\\gt{x,ri)Ma,{k\\)\^ ,  (30) 

where  Mcu  is  the  valence-band  to  conduction-band  dipole  matrix  element 
and  fc||  =  kxX  is  the  wave  vector  along  the  unconfined  direction  of  the  wire. 
When  the  fc||  dependence  of  Mcv  is  neglected,  Eq.  (30)  reduces  to  the  simple 
expression 

\M^\^  =  \MAi^  =  0,r,)\\  (31) 

The  exciton  oscillator  strength  per  unit  length  can  be  written  as  follows 


fls  — 


rriohujgo 


■|M, 


X  2 
cv  • 


(32) 


Here  fiugo  is  the  exciton  ground  state  energy  and  m<,  is  the  free  electron 
mass.  The  oscillator  strength  of  the  excitonic  transition  is  proportional  to  the 
square  of  the  valence-band  to  conduction-band  momentum  matrix  element 
Mcv,  and  to  the  probability  of  finding  the  electron  and  hole  in  the  same  unit 
cell,^®  i.e.,  the  square  of  the  exciton  relative-motion  orbital  wave  function 


gt{x  =  0,7]). 

In  Fig.  10  we  present  the  oscillator  strength  of  the  exciton  transition  per 
unit  length  of  the  wire  as  a  function  of  the  wire  width  for  two  cases:  (i)  with 
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dielectric  confinement;  (ii)  without  dielectric  confinement  .  As  one  can  see, 
the  dielectric  constant  discontinuity  (dielectric  confinement  effect)  causes  a 
20  %  increase  in  the  oscillator  strength  for  the  entire  range  of  the  wire  dimen¬ 
sions.  The  overall  behavior  of  the  oscillator  strength  curve  is  in  an  excellent 
agreement  with  the  data  presented  in  Ref.  [25]  for  the  parabolic  model  ap¬ 
proximation.  The  absolute  value  of  the  oscillator  strength  was  calculated  us- 
ing  the  following  material  parameters:  the  fundamental  band  gap  jE/g  =  1.519 
eV,  the  Kane  matrix  element  Fp  =  23  eV  (related  to  momentum  matrix  ele¬ 
ment),  and  the  density  of  carriers  per  unit  area  =  7.89  x  10^^  /cm^.  For 
the  200  A  wire,  the  oscillator  strength  is  about  0.85/nm  without  dielectric 
confinement  which  is  close  to  the  value  of  0.38/nm  found  in  ref.  [25]  for 
the  same  wire  dimensions.  The  slight  difference  is  due  to  different  values  of 
effective  masses  assumed  in  ref.  [25]. 

Fig.  11  shows  the  oscillator  strength  as  a  function  of  a  magnetic  flux  . 
density.  The  magnetic  field  induced  quenching  of  the  dielectric  confinement 
effect  IS  more  pronounced  in  this  figure  than  in  the  plot  for  the  exciton 
binding  energy.  Since 


/is  oc  J  a{hu})d(jO, 


(33) 


it  appears  that  the  integrated  absorption,  a  quantity  easily  measured  exper¬ 
imentally,  is  quite  sensitive  to  magnetic  field. 

The  matrix  element  for  the  exciton-biexciton  transition  is  given 
by  the  overlap  integral  between  the  exciton  and  biexciton  wave  functions^^ 
which  is  in  our  notation 


=  Ma,  dxj^  dyd^^,  (34) 

where  ^  is  the  exciton  wave  function  given  by  Eq.  (11),  $  is  the  biexciton 
wave  function  given  by  Eq.  (22),  dx  =  dxabdx2b,  dy  =  dy^dy^dy^dy^,  dz  = 
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dzadzbdzidz2,  and  E  is  a  square  region  defined  by  wire  lateral  dimensions. 
In  order  to  capture  only  the  internal  motion  of  the  biexciton,  the  following 
conditions  are  imposed  over  the  coordinates  while  performing  the  integration: 
Xio  =  xi  —  Xa  =  0,  Vi  —  Va  =  0,  zl  —  Za  =  0.  This  condition  keeps  first 
hole  (subscript  “a”)  and  first  electron  (subscript  “1”)  at  the  same  place 
while  allowing  the  other  electron  and  hole  to  move  (here  the  adjective  “first” 
should  be  interpreted  loosely  since  of  course  both  electrons  and  holes  are 
indistinguishable  particles). 

We  have  evaluated  the  integral  in  Eq.  (34)  numerically.  Our  results 
showed  that  the  exciton-biexciton  matrix  element  (or  the  corresponding  os¬ 
cillator  strength)  increases  only  slightly  with  introduction  of  dielectric  con¬ 
finement.  In  order  to  obtain  some  intuition  into  the  origin  of  this  weak  depen¬ 
dence,  we  reproduce  here  the  analytical  result  for  this  quantity  obtained  for 
analogous  but  simpler  trial  wave  functions.  In  Ref.  [27]  the  exciton-biexciton 
wave  function  overlap  was  evaluated  analytically  for  simple  “Gaussian-type” 
orbital  trial  functions,  and  the  exciton  -  biexciton  transition  matrix  element 
was  found  to  be  (in  our  notation) 

=  |M„p[2— V5Fl®|s,(0,T!)pof,  (35) 

O-B 

where  r  is  the  inter-hole  distance  in  the  biexciton  and  D  denotes  the  di¬ 
mensionality  of  the  system.  For  quasi  one-dimensional  structures  {D  —  1) 
and  with  no  magnetic  field  present,  the  oscillator  strength  is  proportional 
to  T  and  inversely  proportional  to  the  exciton  radius  77  (since  gt  is  propor¬ 
tional  to  1/y/^).  Because  both  these  quantities  are  about  equally  affected  by 
the  dielectric  confinement,  they  cancel  each  other’s  effect  so  that  the  overall 
oscillator  strength  is  relatively  insensitive  to  dielectric  confinement. 

Finally,  in  Fig.  12  we  present  the  exciton-biexciton  transition  oscillator 
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strength  as  a  function  of  the  magnetic  flux  density  for  two  wire  widths.  The 
sublinear  behavior  of  the  oscillator  strengths  at  high  values  of  the  magnetic 
flux  density  is  due  to  strong  compression  of  the  biexciton  wave  function  for 
.B  >  8  tesla. 

IV.  CONCLUSION 

In  conclusion,  we  have  theoretically  studied  the  effects  of  dielectric  con¬ 
finement  on  the  binding  energies,  radii  and  oscillator  strengths  of  magneto- 
excitons  and  biexcitons.  Dielectric  confinement  increases  the  binding  energy 
and  oscillator  strength  while  decreasing  the  radius.  A  magnetic  field  tends 
to  reduce  the  effect  of  dielectric  confinement  by  squeezing  the  charged  enti¬ 
ties  (electrons  and  holes)  together  into  an  effective  chargeless  single  particle. 
Although  a  magnetic  field  reduces  the  dielectric  enhancement  of  the  binding 
energy,  it  also  introduces  additional  contribution  to  the  binding  energy  be¬ 
cause  it  squeezes  the  electrons  and  holes  together.  As  a  result,  the  net  binding 
energy  increases  monotonically  with  magnetic  field.  The  resultant  effect  of 
quantum  confinement,  dielectric  confinement  and  magnetostatic  confinement 
may  make  both  the  exciton  and  biexciton  binding  energies  exceed  the  room 
temperature  thermal  energy  in  quantum  wires  of  reasonable  dimension  which 
has  important  ramifications  for  the  experimental  observation  of  these  entities 
in  room  temperature  experiments. 
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V.  APPENDIX 


The  electron  -  image-hole  (hole  -  image-electron)  attractive  terms  for 
the  first  electron  (“1”)  and  the  first  hole  (“a”)  composing  the  biexciton  are 
analogous  to  Eq.  (10)  and  can  be  written  as  follows 


6  >  Cyj  65 


){- 


+ 

+ 

+ 


eyj  eyj  +  et  -I-  yl^  -f  (L^  +  zi  +  ZaY 
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(36) 


+  Via  +  {Lz  -  Zi  -  Za)^ 
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.2 
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\/ ^la  +  (Ly  +  yi  +  Ua)^  + 
_ 1 _ 

+  {Ly  -yi-  VaY  + 


}• 


Here  and  throughout  the  Appendix  we  have  used  the  same  notation  for  the 
coordinates  as  in  the  text.  The  remaining  attraction  terms 
and  UattT^^  can  be  obtained  by  formal  substitution  of  appropriate  indices  in 
Eq.  (36). 

The  self-energy  terms  are  analogous  to  Eqs.  (8),  (9)  and  for  the  first  hole 
are  given  by 


rr(a) 

^self 


I  1  C  .  £tv  )  f  ^  ^ 

2 ett)  €10  +  66  |2za  —  Til  \2za-\-L2;\ 
1  1 


(37) 


'*'|2y. +  L,r 

The  rest  of  the  self-energy  terms,  ufjlf,  and  can  be  obtained  by 
formal  substitution  of  appropriate  indices  in  Eq.  (37). 

The  only  new  terms  which  appear  in  the  biexciton  problem,  but  do  not 
appear  in  the  exciton  problem,  are  the  repulsive  terms  between  electron  - 
“another  electron  image”  (hole  -  “another  hole  image”).  These  terms  have 
the  same  structure  as  the  attractive  terms  but  with  the  opposite  sign.  The 
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repulsive  term  for  the  electrons  can  be  written  as 


^  repuls 


+ 

+ 

+ 


1 


\Jxi2  +  1/12  +  {Lz  -Z\-  22)^ 
_ 1 _ 

y/ ^12  +  {Ly  +  yi+  2/2)^  +  2:12 


^ XI2  +  {Ly  -yi-  2/2)^  +  2^12 


(38) 


Again,  the  rest  of  the  repulsive  terms  can  be  obtained  by  formal 

substitution  of  appropriate  indices  in  Eq.  (38). 

The  group  of  attractive  terms  contains  16  x  2  separate  interactions,  the 
group  of  self-energy  terms  has  16  interactions  and  the  group  of  the  repulsive 
terms  has  8x2  separate  interaction.  Combining  these  terms  together,  we  have 
64  terms  accounting  for  the  Coulomb  interaction  between  all  charged  entities 
composing  the  biexciton  and  their  images  in  the  four  plane  boundaries  of  the 
wire. 
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FIGURE  CAPTIONS 


[1] .  (a)  The  geometry  of  the  wire,  (b)  Schematic  structure  of  the  wire  along 
y  direction  and  illustration  of  the  image-charge  method.  The  original  charge 
is  indicated  by  a  double  circle  and  the  induced  image  charges  are  given  by 
single  circles. 

[2] .  The  magneto-electric  subband  bottom  energies  and  as  func¬ 
tions  of  the  wire  width  for  two  different  values  of  a  magnetic  flux  density.  The 
wire  thickness  is  200  A.  The  two  upper  curves  correspond  to  the  electrons 
and  the  two  lower  to  the  holes.  The  curves  delineated  by  “o”  are  plotted  for 
zero  magnetic  flux  density  and  those  delineated  by  “-t-”  are  plotted  for  a  10 
tesla  magnetic  field.  The  fundamental  band  gap  energy  Eq  is  assumed  to  be 
zero  for  clarity. 

[3] .  Exciton  binding  energy  in  a  GaAs  quantum  wire  as  a  function  of  the  wire 
width  Ly  for  three  different  cases;  (i)  without  dielectric  confinement  (dashed 
line);  (ii)  with  only  first  order  correction  due  to  the  dielectric  confinement 
(dash-dotted  line);  and  (iii)  with  complete  dielectric  confinement  (solid  line). 
In  the  last  case,  up  to  fourth  order  correction  was  sufficient  for  convergence. 
The  thickness  of  the  wire  along  the  z-direction  is  200A. 

[4] .  Exciton  longitudinal  radii  versus  wire  width  with  and  without  dielectric 
confinement.  When  there  is  no  magnetic  field  present,  the  exciton  radius 
increases  monotonically  with  increasing  wire  width.  The  thickness  along  the 
z-direction  is  200A. 

[5] .  Exciton  longitudinal  radii  as  a  function  of  a  magnetic  flux  density  with 
(“-f”)  and  without  (“o”)  dielectric  confinement.  The  thickness  along  the 
z-direction  is  200A. 
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[6] .  Exciton  binding  energy  versus  magnetic  field.  All  parameters  are  the 
same  as  in  Fig.  5.  The  inset  shows  exciton  binding  energy  as  a  function  of 
the  dielectric  misfit. 

[7] .  Exciton  radius  as  a  function  of  the  dielectric  misfit  for  two  values  of  the 
magnetic  field:  (i)  B=0  tesla  (“o”)  and  (ii)  B=10  tesla  (“+”)•  The  thickness 
along  the  ^-direction  is  200A  and  the  width  is  500A. 

[8] .  Exciton  binding  energy  as  a  function  of  the  dielectric  misfit  for  two 
values  of  the  magnetic  field:  (i)  B=0  tesla  (“o”);  (ii)  B=10  tesla  (“+”).  The 
thickness  along  the  ^-direction  is  200A  while  the  width  is  500A. 

[9] .  (a)  Biexciton  binding  energy  in  a  GaAs  quantum  wire  as  a  function  of 
magnetic  flux  density,  (b)  Biexciton  binding  energy  in  a  GaAs  quantum  wire 
as  a  function  of  wire  width.  No  magnetic  field  is  present.  The  thickness 
along  the  z-direction  is  200A  and  the  width  is  500A  for  both  cases. 

[10.  Exciton  oscillator  strength  as  a  function  of  wire  width  with  (“+”)  and 
without  (“o”)  dielectric  confinement.  The  thickness  along  the  z-direction  is 

200A. 

[11] .  Exciton  oscillator  strength  as  a  function  of  a  magnetic  flux  density  with 
(“-h”)  and  without  (“o”)  dielectric  confinement.  The  thickness  along  the 
z-direction  is  200 A  and  the  width  is  500 A. 

[12] .  The  exciton-biexciton  transition  oscillator  strength  as  a  function  of  a 
magnetic  flux  density.  For  the  curve  delineated  by  “o”  the  thickness  along 
the  z-direction  is  200A  and  the  width  is  500A.  For  the  curve  delineated  by 
“+”  the  thickness  along  the  z-direction  is  200A  and  the  width  is  700A. 
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We  have  calculated  the  momentum  and  energy  relaxation  rates  of  hot 
electrons  in  a  quantum  wire  subjected  to  a  magnetic  field.  These  rates 
have  been  determined  for  both  forward  (where  the  electron  velocity  does 
not  change  sign)  and  backward  (where  the  electron  velocity  changes  sign) 
scattering.  A  magnetic  field  suppresses  both  momentum  and  energy  relax¬ 
ation  rates  associated  with  non-polar  acoustic  and  surface  optical  phonon 
scattering.  For  polar  acoustic  and  polar  optical  phonon  scattering,  it 
suppresses  the  momentum  relaxation  rate,  but  increases  the  energy  re¬ 
laxation  rate.  The  suppression  of  momentum  relaxation  rate  leads  to  an 
increase  in  the  drift  velocity  and  mobility  of  electrons,  which  is  beneficial 
for  electronic  device  applications. 


1.  Introduction 

It  is  now  well-established  that  a  magnetic  field  suppresses  back-scattering  of 
electrons  in  quantum  wires.  The  field  spatially  separates  oppositely  traveling 
electron  states  (  edge-states  )  thereby  reducing  the  overlap  between  their  wave 
functions  and  the  matrix  element  for  backscattering  [1-4].  While  this  has  many 
important  consequences,  perhaps  the  most  noteworthy  is  the  celebrated  integral 
quantum  Hall  effect  whose  origin  may  lie  in  this  phenomenon  [5]. 

In  this  paper,  ,we  have  calculated  the  momentum  and  energy  relaxation  rates 
of  electrons  in  a  quantum  wire  subjected  to  a  magnetic  field.  Since  the  mo¬ 
mentum  relaxation  rate  is  primarily  determined  by  backscattering  (large  change 
in  momentum)  as  opposed  to  forward  scattering  (small  change  in  momentum), 
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We  have  theoretically  studied  giant  dipoles  associated  with  transitions 
between  magneto-electric  subbands  in  a  quantum  wire  subjected  to  a  trans¬ 
verse  magnetic  field.  The  strengths  of  these  dipoles  and  their  resonant  fre¬ 
quencies  can  be  varied  with  the  magnetic  field  which  then  allows  one  to  tune 
the  emission  wavelength  of  these  transitions.  The  large  magnitude  of  the 
dipole  moments  also  leads  to  a  strong  second-harmonic  component  of  the  di¬ 
electric  susceptibility  which  can  be  utilized  for  non-linear  optical  applications 
such  as  second-harmonic  generation,  limiting,  mixing,  optical  switching,  etc. 
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I.  Introduction 


Direct  intraband  transitions  between  the  quantized  states  (subbands)  of 
the  conduction  band  in  a  quantum  well  is  a  well-researched  topic. ^  It  has 
been  shown  both  experimentally  and  theoretically  that  such  transitions  have 
very  large  dipole  moments  and  narrow  bandwidths.  Strong  infrared  absorp¬ 
tion  associated  with  transitions  between  the  lowest  two  electronic  subbands 
in  a  GaAs  quantum  well  was  observed  long  ago  by  a  number  of  experimental 
groups^.  Recently,  population  inversion  between  the  second  and  third  sub¬ 
bands  of  a  quantum  well  has  been  established  unambiguously  leading  to  the 
demonstration  of  the  celebrated  quantum  cascade  laser^.  The  energy  sepa¬ 
ration  between  the  subbands  in  a  quantum  well  or  wire  can  be  varied  by  an 
external  magnetic  field  which  then  allows  one  to  realize  a  continuously  tun¬ 
able  laser  or  light-emitting-device.  Moreover,  the  field  can  induce  forbidden 
transitions  that  make  additional  frequency  ranges  accessible,  thus  permitting 
flexible  device  design. 

Another  potential  use  of  magnetic  field  biasing  of  quantum  wells  or  wires 
is  in  non-linear  optics.  Non-linear  optical  properties  stem  from  higher  or¬ 
der  dielectric  susceptibilties.  Specifically,  the  second-order  susceptibility 
is  responsible  for  such  phenomena  as  mixing  and  second-harmonic  genera¬ 
tion.  It  is  well-known  that  even-order  susceptibilities  vanish  in  structures 
with  inversion  symmetry.  Consequently,  finite  second-order  susceptibilities 
can  be  obtained  in  such  structures  only  if  the  inversion  symmetry  of  the 
conduction-band  potential  is  broken  either  by  an  external  electric  field  or  by 
the  intentional  growth  of  an  aisymmetric  well.  Obviously,  the  former  is  the 
preferred  method  since  an  electric  field  can  be  continuously  varied  which  al¬ 
lows  one  to  tune  the  degree  of  symmetry-breaking  and  the  magnitude  of 
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However,  this  method  has  a  practical  shortcoming.  An  electric  field  tilts  the 
potential  barriers  of  the  well  thereby  allowing  carriers  to  escape  by  tunneling 
or  thermionic  emission.  This  is  especially  serious  in  GaAs/AlGaAs  systems 
where  the  barrier  height  is  relatively  small.  Indeed,  it  has  been  pointed  out 
that  the  electronic  states  in  a  quantum  confined  system  biased  by  a  trans¬ 
verse  electric  field  are  never  true  bound  states  since  the  particles  can  always 
lower  their  energy  by  escaping  from  the  welH.  Therefore,  these  states  have  a 
finite  lifetime.  An  electric  field  reduces  the  lifetime  drastically  and  broadens 
the  transitions. 

To  overcome  this  shortcoming,  one  can  adopt  magnetostatic  biasing.  A 
magnetic  field  can  break  inversion  symmetry  without  tilting  potential  barri¬ 
ers  and  promoting  carrier  escape.  A  transverse  magnetic  field  applied  to  a 
quantum  wire  exerts  a  Lorentz  force  on  an  electron  moving  along  the  wire. 
As  a  result,  its  wave  function  (in  any  subband)  will  be  skewed  towards  one 
edge  of  the  wire.  This  skewing  does  not  tilt  potential  barriers  to  first  order 
(the  barriers  may  tilt  slightly  because  of  a  second-order  effect  associated  with 
space-charges  and  the  self-consistent  (Hall)  electric  field).  However,  it  effec¬ 
tively  breaks  inversion  symmetry  since  it  causes  a  net  charge  to  accumulate 
at  either  edge  of  the  wire  (the  charges  at  the  two  edges  have  opposite  signs  as 
in  the  classical  Hall  effect) .  This  leads  to  a  non- vanishing  even-order  suscep¬ 
tibility  in  a  symmetric  structure.  The  skewing  has  another  subtle  effect.  The 
degree  to  which  the  wave  function  is  skewed  is  different  in  different  suhbands 
since  an  electron  has  different  kinetic  energies  and  hence  experiences  differ¬ 
ent  Lorentz  forces  in  different  subbands.  As  a  result,  transitions  between 
subbands  whose  wave  functions  have  the  same  parity  -  which  are  forbidden 
without  a  magnetic  field  -  are  now  allowed  since  the  parities  are  altered  by 
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different  amounts  in  different  subbands  by  the  different  degrees  of  skewing. 
We  should  point  out  that  this  effect  has  striking  similarity  with  the  quan¬ 
tum  confined  Lorentz  effect  (QCLE)  previously  examined  by  Balandin  and 
Bandyopadhyay  ®  in  the  context  of  interband  transitions  between  conduction 
and  valence  band  states  rather  than  intraband  transitions  in  the  conduction 
band.  In  that  case,  the  different  amount  of  skewing  of  the  electron  and  hole 
states  quenches  photoemission. 

This  paper  is  organized  as  follows.  In  the  next  section,  we  describe  the 
theoretical  formulation,  followed  by  results.  Finally,  in  section  IV,  we  present 
the  conclusions. 

II.  Theory 

We  consider  a  quantum  wire  as  shown  in  Fig.  1  with  a  magnetic  field 
applied  along  the  z  direction.  The  thickness  along  the  2:  direction  is  so  small 
(and  consequently  the  subband  separation  in  energy  is  so  large)  that  for 
the  range  of  photon  energies  considered,  an  electron  cannot  be  excited  (by 
real  transition)  into  a  subband  which  has  more  than  two  nodes  along  the  z- 
direction.  Such  a  transition  will  not  be  accessible  in  energy.  This  restriction, 
coupled  with  the  fact  that  a  magnetic  field  does  not  affect  the  z-component 
of  the  electron  wave  function,  allows  us  to  drop  the  z-component  from  further 
consideration.  The  width  of  the  wire  along  the  y-direction  is  however  large 
enough  that  subbands  with  more  than  two  nodes  along  the  y-direction  are 
accessible  in  energy. 

In  the  framework  of  the  envelope  function  approximation  (EFA),  an  elec¬ 
tron  wave  function  can  be  written  as  the  product  of  a  Bloch  wave  function, 
periodic  with  the  atomic  lattice  spacing,  and  an  envelope  wave  function,  de¬ 
scribing  the  nonperiodic  behavior.  Consequently,  the  wave  function  of  an 
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electron  for  a  given  wave  vector  k  along  the  a:-direction,  in  the  nth  magne¬ 
toelectric  subband,  at  a  magnetic  field  B  can  be  written  as 

^n{x,  y,  k,  B,  t)  =  y,  k,  B)un{x,  y,  (1) 

where  y,  B)  is  an  envelope  function,  n„(x,  y,  k)  is  a  Bloch  function  of 
a  conduction  band  and  En{k,  B)  is  the  dispersion  relation  of  the  nth  magne¬ 
toelectric  subband  at  a  flux  density  B.  The  Bloch  wave  functions  are  assumed 
to  be  s  states  which  is  the  usual  case  for  semiconductors  where  J  =  1/2  for 
the  conduction  band. 

The  envelope  function  can  be  further  decomposed  into  a  plane  wave  along 
the  unconfined  x  direction  and  confined  component  along  y  direction 

<6^{x,y,k,B)  =  Xn(y,k,By‘“  (2) 

Using  the  electric  dipole  approximation,  we  can  write  the  matrix  element 
of  photoinduced  inter-subband  transitions  within  the  conduction  band  as® 

dj^i{k,  B)  =  eJ  Xf{y,  k,  B)t)  •  fx.(y,  k,  B)dr  J  u){x,  y,  k)ui{x,  y,  k)d^l  (3) 

where  dCl  is  a  volume  element,  rj  is  the  unit  vector  along  the  direction  of 
the  incident  photon  polarization,  f  =  xa^  -b  yoy  is  the  two-dimensional  ra¬ 
dius  vector,  and  subscripts  i,  f  stand  for  initial  and  final  states  respectively. 
The  exponential  term  of  Equation  (2)  is  not  present  in  Eq.  (3)  since  for 
photoinduced  transitions  {kf  =  A;,),  the  product  of  the  exponential  function 
and  its  complex  conjugate  is  exactly  unity.  The  volume  overlap  of  the  Bloch 
functions  is  also  unity  for  s-states  with  the  same  wave  vector.  Now,  if  we 
assume  that  the  incident  light  is  polarized  along  the  y-direction  so  that  t}  = 
dy,  the  above  equation  simplifies  to 
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rWl2 

=  e  <  Xf\y\Xi  >=  e  /  yxAy,k,B)xi{y,h,B)dy,  (4) 

J-WI2 


where  W  is  the  width  of  the  quantum  wire  along  the  y-direction. 

One  should  note  here,  that  if  there  is  no  magnetic  (or  electric)  field  ap¬ 
plied,  the  envelope  functions  Xi  just  particle-in-box  states  and  the  dipole 
moment  in  Eq.  (4)  is  non  zero  only  for  the  transitions  between  subband 
states  of  opposite  parity.  For  a  symmetric  square  potential  well,  these  dipole 
elements  (between  any  two  states  n  and  m)  are  independent  of  the  wave  vec¬ 
tor  k  and  can  be  found  analytically^  by  evaluating  the  integral  in  Equation 

(4). 


di,i  =  e  <  Xf\y\Xi  > 


eW-T-. — - if  n  and  m  have  opposite  parity 

0,  otherwise  (5) 


However,  when  a  magnetic  field  is  applied,  the  skewing  of  the  wave  functions 
changes  the  integral  in  Equation  (4)  and  alters  the  selection  rules.  Generally, 
the  skewing  causes  three  effects.  First,  it  makes  the  dipole  moment  depend 
on  the  wave  vector  k  (since  the  degree  of  skewing  depends  on  k).  Second,  it 
reduces  the  dipole  moment  for  transitions  between  states  of  opposite  parity 
(since  the  integral  in  Equation  (4)  decreases),  and  third,  it  allows  forbidden 
transitions  between  states  of  the  same  parity  (since  the  integral  in  Equation 
(4)  no  longer  vanishes  for  states  of  the  same  parity). 

It  is  clear  from  Equation  (4)  that  to  calculate  the  dipole  moments  in  the 
presence  of  a  magnetic  field,  all  we  need  to  compute  are  the  wave  functions 
X/,»(y>  k,  B)  at  a  given  magnetic  field  B,  for  given  magnetoelectric  subbands 
/  and  i,  and  for  a  given  wave  vector  k.  This  is  achieved  via  a  numerical  (finite 
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difference)  solution  of  the  Schrodinger  equation  following  the  prescription  of 
ref.  [7].  Once  this  is  done,  we  can  calculate  the  dipole  moment  in  Equation 
(4)  for  any  chosen  intersubband  transition,  at  any  chosen  magnetic  field  and 
for  any  chosen  wave  vector. 

In  the  limit  of  high  magnetic  fields,  when  the  magnetic  length  l(= 
eB)  «W  one  can  again  obtain  an  analytical  expression  for  the  dipole 
moment  d/,;.  In  this  case,  the  magnetostatic  confinement  predominates  over 
electrostatic  confinement  and  the  envelope  functions  Xn{yik,B)  can  be  ap¬ 
proximated  by  harmonic-oscillator  wave  functions: 

X„{y,k,B)  =  =  KH^(a,y  -  (6) 

where  iV„  =  is  a  normalization  constant,  Hn{oc,y)  is  the  nth 

Hermite  polynomial,  yk  =  hk/eB,  and 

/^_  1 
V  “  /■ 

In  order  to  evaluate  the  integral  in  Eq.  (4)  analytically,  we  extend  the  lim¬ 
its  of  integration  to  infinity  assuming  that  the  wave  function  tail  is  negligible 
at  the  boundaries  of  the  wire  (i.e.  at  y  =  ±17/2).  This  is  a  very  reason¬ 
able  assumption  in  a  high  confining  magnetic  field.  The  resulting  analytical 
expression  for  the  dipole  moment  is 

df,i{B)  =  e<xAy\Xi>  =  ifm  =  n  +  l  (8) 

=  if  m  =  n-l 

=  0,  otherwise. 

The  physical  significance  of  the  two  analytical  limits  -  5  — >  0  and  5  — >  oo 
-  is  obvious.  At  zero  field,  the  dipole  is  determined  by  the  width  of  the 
wire  d/,,-  ~  eW,  and  at  the  high  field  limit  it  is  determined  by  the  magnetic 
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length  df^i  ~  el.  This  is  what  one  would  expect  intuitively.  At  zero  field, 
the  dipole  is  confined  electrostatically  with  the  wire  width  being  a  measure 
of  this  confinement  while  at  high  magnetic  field,  the  dipole  is  confined  mag- 
netostatically  and  the  magnetic  length  is  the  corresponding  measure  of  this 
confinement, 

III.  Results 

A.  Intraband  dipoles 

We  now  present  results  of  our  calculations.  The  physical  parameters 
used  for  the  numerical  calculations  correspond  to  a  GaAs  quantum  wire  with 
relative  dielectric  constant  =  12.9,  and  effective  masses  me  =  0.067mo  and 
TTih  =  O.Stoo  where  mo  is  the  free  electron  mass. 

In  Fig.  2,  we  show  the  dependence  of  the  dipole  moment  df,i{k,B)  for 
three  transitions  (el-e2,  e2-e3,  and  el-e3)  on  the  wave  vector  k  when  a  mag¬ 
netic  field  of  1  tesla  is  applied  (following  usual  practice,  the  transitions  are 
numbered  by  the  subband  indices).  The  dipoles  corresponding  to  transitions 
between  states  of  opposite  parity  (el-e2  and  e2-e3)  have  maxima  at  A:  =  0 
and  then  decrease  with  increasing  wave  vector.  This  can  be  easily  understood 
as  follows.  At  zero  wave  vector  (no  translational  velocity)  these  states  do  not 
experience  any  Lorentz  force  and  hence  the  wave  functions  are  not  skewed. 
As  the  wave  vector  k  increases,  the  translational  velocity  and  the  Lorentz 
force  experienced  increase.  Consequently,  the  envelope  wave  functions  are 
skewed  more  and  more  and  the  dipole  moment  decreases.  Real  transitions 
between  states  of  the  same  parity  are  forbidden  at  zero  magnetic  field,  but 
at  a  finite  magnetic  field,  they  are  forbidden  only  at  A;=0  when  there  is  no 
translational  velocity  and  no  Lorentz  force  to  skew  the  wave  functions.  With 
increasing  k,  the  wave  functions  are  increasingly  skewed  and  the  dipole  mo- 
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ment  of  forbidden  transitions  increases.  In  our  chosen  prototype  wire,  dei-es 
reaches  a  maximum  of  28  e  —  A  at  ^  =  0.0051  A  and  then  decreases  gradually 
ultimately  reaching  zero.  This  intriguing  non-monotonic  dependence  on  k  is 
explained  later  on.  However,  at  this  point,  it  is  interesting  to  note  that  a 
fairly  large  forbidden  dipole  moment  of  rsj  30  e-A  can  be  achieved  in  realistic 
structures  at  a  moderate  magnetic  field  of  1  tesla. 

Fig.  3  presents  the  dipole  moments  for  the  same  transitions  cis  a  func¬ 
tion  of  magnetic  flux  density.  The  propagation  wave  vector  k  is  chosen  to 
be  0.01/A.  At  zero  magnetic  field,  a  non- vanishing  dipole  matrix  element 
occurs  only  for  transitions  between  states  of  opposite  parity  (el-e2,  e2-e3) 
as  expected  from  Equation  (5).  This  equation  also  allows  us  to  estimate  the 
strengths  of  these  zero-field  dipoles  to  be  180  e  —  A  for  el-e2  and  195  e  —  A 
for  e2-e3  transitions.  As  we  can  see  from  Fig.  3,  these  values  are  in  excel¬ 
lent  agreement  with  our  numerical  results.  From  the  analytical  expression  in 
equation  (8),  we  can  estimate  the  strength  of  the  el-e2  dipole  to  be  66  e  —  A 
at  a  magnetic  flux  density  of  15  tesla.  This  number  also  agrees  with  our 
numerical  result.  The  el-e3  dipole  vanishes  at  both  zero  field  (because  of  the 
spatial  symmetry  of  the  particle-in-a-box  states)  and  at  high  fields  because 
of  the  symmetry  of  the  Landau  states  or  Hermite  polynomials.  This  behav¬ 
ior  is  consistent  with  Equations  (5)  and  (8).  Only  at  intermediate  fields, 
when  the  wave  functions  of  the  subbands  are  a  hybrid  between  particle-in- 
a-box  states  and  Hermite  polynomials  (and  thus  “non-symmetric”  in  space), 
is  this  transition  allowed.  This  immediately  tells  us  that  dgz-ei  must  have 
a  non-monotonic  dependence  on  the  magnetic  flux  density  B  and  indeed  it 
does. 

A.l  Non-monotonic  behavior  of  dipole  moment  des-ei* 
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Let  us  now  examine  the  non-monotonic  behavior  of  rfes-ei  more  closely. 
This  transition  is  forbidden  at  zero  field  since  the  wave  functions  of  the 
first  and  third  subband  have  the  same  parity.  At  low  magnetic  fields,  the 
parities  are  altered  by  the  skewing  of  the  wavefunctions  and  consequently 
de3_ei  is  no  longer  zero  but  increases  with  the  magnetic  field.  It  reaches 
a  maximum  of  about  30  e-A  and  then  decreases.  This  later  decreases  is 
related  to  the  following  effect.  For  a  fixed  wave  vector  k,  a  sufficient  increase 
in  the  flux  density  B  forces  the  traversing  states  (“skipping  orbits”  or  “edge 
states”)  to  condense  into  closed  cyclotron  orbits  (Landau  levels)  which  are 
no  longer  skewed  by  the  magnetic  field  to  the  wire  edge  since  they  have  no 
translational  velocity  and  hence  no  Lorentz  force.  While  edge  states  have 
a  skewed  wave  function  which  is  not  symmetric  in  space,  cyclotron  orbits 
have  a  wave  function  that  is  symmetric  about  the  orbit  center  yk.  Note  that 
yk  depends  only  on  k  and  B.  Therefore,  at  a  fixed  k,  the  wave  functions 
of  the  first  and  third  Landau  levels  are  symmetric  about  a  common  center. 
Whenever  this  kind  of  symmetry  holds,  de3_ei  vanishes.  Therefore,  the  dipole 
moment  dgs-ei  decreases  gradually  to  zero  at  high  magnetic  field  with  the 
onset  of  Landau  condensation. 

The  same  physics  can  be  elucidated  in  a  different  way  by  considering  the 
energy  versus  wave  vector  relation  is  Fig.  4(a)  and  4(b)  which  show  the 
dispersion  of  the  first  and  third  magneto-electric  subbands  respectively. 

At  5  =  0,  velocity  (slope  of  the  curves)  aX  k  =  0.01  /A  are  non-zero  for 
both  the  el  and  e3  subbands.  However,  the  Lorentz  force  is  zero  because 
B  =  0  and  hence  dei-ez  =  0.  At  B  =  5  tesla,  the  group  velocities  for  the 
two  subbands  are  still  non-zero  and  the  Lorentz  force  is  finite  resulting  in 
skewing  of  wave  functions  and  a  non-vanishing  value  of  dgx-ez-  At  jB  =  10 
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tesla,  the  group  velocities  at  fc=0.01  /A  are  zero  in  both  subbands  indicating 
that  the  corresponding  states  have  undergone  Landau  condensation.  In  this 
case,  the  Lorentz  force  (for  skewing)  is  again  zero  and  the  dipole  moment 
dei-e3  vanishes  once  more.  The  crucial  point  to  note  is  that  the  Lorentz 
force  ev  X  B  can  vanish  in  two  different  ways:  (i)  B  =  0,  and  (ii)  u  =  0. 
These  two  conditions  are  met  at  zero  and  very  high  magnetic  fields.  As 
a  result,  the  dipole  moment  dei-es  exhibits  a  non- monotonic  behavior  in 
magnetic  field.  One  can  ask  as  to  why  the  same  physics  does  not  cause 
non-monotonicity  in  the  el-e2  and  e2-e3  curves.  It  is  not  clear  apriori  that 
non-monotonicity  cannot  occur  (indeed  there  are  regions  of  inflexion  in  the 
two  curves).  However,  the  point  to  note  is  that  Landau  condensation  causes 
recovery  of  the  wave  function  symmetry  (or  anti-symmetry),  but  does  not 
restore  the  original  zero-field  wavefunctions.  This  is  shown  in  Fig.  5  where 
we  show  the  wave  functions  in  el  subband  at  0  and  10  tesla.  Both  wave 
functions  are  “symmetric”  in  space,  but  they  are  otherwise  vastly  different 
since  the  magnetostatic  confinement  squeezes  the  wave  functions  binding 
them  in  cyclotron  orbits. 

The  non-monotonicity  in  the  wave  vector  dependence  of  dei_e3  in  Fig. 
2  has  a  similar  origin.  As  k  is  increased,  the  relative  skewing  between  the 
wave  functions  in  el  and  e3  subbands  change  non-monotonically  causing  the 
non-monotonicity  in  this  figure. 

The  process  described  above  is  illustrated  in  Fig.  6  (a)  -  (c),  where  we 
present  wave  functions  of  two  electronic  states  (el  and  e3)  for  three  values  of 
magnetic  flux  density.  At  zero  magnetic  field  the  wave  functions  are  symmet¬ 
ric  about  the  center  of  the  wire  and  dipole  transition  dez-ei  is  forbidden  (case 
“a”).  At  low  magnetic  field  the  wave  functions  are  skewed  to  the  edge  of  the 
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wire  (“edge  states”  -  case  “b”)  and  the  spatial  symmetry  is  broken  for  both 
states*  Dipole  transition  dgs-ei  is  now  allowed.  It  is  important  to  note  here 
that  the  symmetry  breaking  skewing  of  the  wave  functions  is  caused  by  the 
simultaneous  presence  of  a  magnetic  field  and  the  electrostatic  potential  bar¬ 
riers  at  the  edges  of  the  quantum  wire.  At  higher  magnetic  fields,  when  the 
magnetic  length  is  smaller  than  the  wire  width,  the  electrons  do  not  “feel”  the 
potential  barriers  at  the  edges  of  the  wire  as  they  undergo  complete  Landau 
condensation  and  execute  cylotron  motion  with  radius  much  smaller  than 
the  width  of  the  wire.  In  this  case,  the  wave  function  symmetry  is  essentially 
restored  (case  “c”)  although  the  wave  functions  are  now  symmetric  about  a 
point  that  is  not  at  the  center  of  the  wire.  Nonetheless,  what  is  important  is 
that  both  wave  functions  are  symmetric  about  the  same  point.  Consequently, 
the  c?e3-ei  transition  vanishes.  Simultaneous  presence  of  both  electrostatic 
confinement  and  magnetostatic  confinement  is  therefore  necessary  for  wave 
function  skewing,  formation  of  edge  states  and  the  observation  of  forbidden 
transitions. 

B.  Second  Harmonic  Generation 

It  is  weU  known  that  in  systems  with  inversion  symmetry  there  can  be  no 
second  order  nonlinearity.®  However,  in  systems  without  inversion  symmetry, 
the  lowest  order  optical  nonlinearity  is  of  the  second  order  and  is  expressed 
by 

pT2)(fc,a;)  =  x^'^K^',^i,^2)Ei{ki,u;i)E2{k2,uj2),  (9) 

where  P  is  the  polarization  caused  by  two  electric  fields  Ei  and  E2  that  are 
associated  with  the  electromagnetic  fields  of  either  two  frequency  components 
of  the  same  light  beam  or  two  different  coherent  beams  with  frequencies  a?,- 
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and  wave  vectors  fc,-.  The  frequencies  and  wave  vectors  obey  the  energy  and 
momentum  conservation  laws 


hu!  =  (10) 

i 

hk  =  ^  iiUki, 

t 

It  is  obvious  that  the  third-ranked  tensor  will  vanish  in  any  structure 
with  inversion  symmtery.  A  quantum  confined  structure  may  lack  inver¬ 
sion  symmetry  for  two  main  reasons,  (i)  the  semiconductor  material  by  its 
intrinsic  chemical  and  crystalline  structure  may  lack  inversion  symmetry,® 
and  this  is  the  case  in  most  IITV,  ITVI,  and  I- VII  compounds  along  certain 
crystallographic  directions,  (ii)  the  quantum  confining  potential  well  may  be 
asymmetric  (e.  g.  triangular  potential  well,  asymmetric  double  square  well 
potential,  etc.).  In  the  first  case,  the  asymmetry  is  related  to  the  intracell 
charge  asymmetry  and  is  not  affected  by  the  the  confinement  since  the  latter 
extends  over  several  unit  cells.  In  the  second  case,  the  asymmetry  is  arti¬ 
ficially  imposed  and  therefore  can  be  engineered.  It  clearly  depends  on  the 
confining  potential  and  hence  an  applied  electric  field  can  alter  the  potential 
and  change  the  degree  of  symmetry-breaking. 

In  the  present  work  we  restrict  ourselves  to  the  second  case  and  do  not 
consider  intrinsic  second  order  nonlinearities  of  GaAs  which  are  actually  quite 
large  (the  nonlinear  susceptibility  of  bulk  GaAs  is  10“^°  m/V.^®  ) 

As  mentioned  before,  we  avoid  an  electric  field  since  it  promotes  carrier 
escape  and  consider  a  magnetic  field  instead.  Although  a  magnetic  field  does 
not  affect  the  potential  to  first  order,  the  simultaneous  action  of  symmetric 
electrostatic  potential  and  an  external  magnetic  field  may  lead  to  the  uneven 
charge  distribution  along  the  width  {y  —  axis)  of  the  wire  caused  by  different 
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degrees  of  skewing  of  the  wave  functions.  Because  of  this  reason,  it  is  possible 
to  break  the  inversion  symmetry  in  a  symmetric  quantum  well  or  wire  with 
a  magnetic  field  alone.  This  approach  is  superior  to  applying  a  transverse 
electric  fields  since  the  latter  will  tilt  the  confining  potential  wells  thereby 
promoting  carrier  escape  from  the  well  by  either  tunneling  or  thermionic 


emission. 


The  large  magnitude  of  the  dipole  moments  associated  with  otherwise 
forbidden  transitions  between  subbands  of  the  same  parity  and  their  sensi¬ 
tivity  to  the  biasing  magnetic  field  open  up  a  possibility  of  second  harmonic 
generation  (SHG)  that  can  be  controlled  by  the  magnetic  field.  In  order  to 
evaluate  the  magnitude  and  dependences  of  SHG  on  the  biasing  field  and  wire 
geometry,  we  calculate  the  second  order  susceptibility  using  the  formula^^ 


abc 


dah^bAa 


{Qba  -U>1-  U!2){^ca  “  ^2) 


1,  (11) 


where  N  is  concentration  (number  density)  of  conduction  electrons,  hQap  = 
hUap{B ,W)  is  the  energy  spacing  between  a,P  subbands  which  depend  on 
the  applied  magnetic  field  and  wire  width,  dmn  =  dmn{B,  W)  is  a  dipole  ele¬ 
ment  calculated  using  Eq.  (4),  and  is  defined  to  be  The  total 

symmetrisation  operation  St  indicates  that  the  expression  which  follows  it  is 
to  be  summed  over  all  six  permutations  of  the  pairs  (ju,  —uJa),  (/?j‘^2)- 

Since  Sj  involves  a  summation  over  all  possible  permutations,  it  is  clear  that 
is  invariant  under  any  of  them.  For  simplicity,  the  Fermi 
distribution  po{o)  was  assumed  to  be  unity. 

Eq.  (11)  is  an  approximation  which  applies  only  under  the  condition  that 
all  of  the  optical  frequencies  involved  (operational  frequencies  are 

removed  far  enough  from  the  subband  transition  frequencies.  It  means  that 
the  medium  is  assumed  to  be  transparent  and  loss-free  at  all  the  relevant 
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optical  frequencies.  This  assumption  can  be  relaxed  by  the  introduction  of 
transition  damping  factors  into  the  expression  in  Eq.  (11).  In  our  consid¬ 
eration  we  are  mainly  interested  in  the  effects  of  an  applied  magnetic  field 
on  the  second  order  susceptibility.  Since  these  effects  manifest  themselves  in 
Eq.  (11)  primarily  via  the  magnetic  field  dependence  of  the  dipole  elements 
dmn  =  dmn{B)i  we  did  not  include  any  damping  constants  and  associated 
finite  linewidths  of  the  electronic  states.  One  should  also  note  here  that  the 
Eq.  (11)  is  strictly  correct  only  for  dilute  media.  In  this  case,  one  can  write 
^(2)  =  with  being  the  2nd-order  nonlinear  polarization.  The  above 

expression  is  valid  only  under  moderate  excitation. 

In  Fig.  7,  we  present  normalized  values  of  as  a  function  of  magnetic 
field  for  three  different  wire  widths  and  a  fixed  value  of  the  wave  vector  k 
(fixed  excitation  frequency).  The  operational  frequencies  ui  —  are  chosen 
for  a  CO2  laser.  For  wide  ranges  of  magnetic  flux  densities  {B  <  20  tesla) 
and  wire  widths  (lOOA  <  W  <  lOOOA),  these  frequencies  are  removed  far 
enough  from  the  subband  transition  frequencies  As  long  as  the 

latter  is  true,  the  x^^^  dependence  on  magnetic  field  is  governed  mainly  by 
dipole  elements  dmn-  Consequently,  the  x^^^  curve  for  W  =  lOOOA  peaks 
at  the  same  value  of  a  magnetic  flux  density  (B=0.3  tesla)  as  the  el-e3 
dipole  curve  of  Fig.  3.  The  magnetic  flux  density  at  which  x^^^  reaches 
its  maximum  increases  with  decreasing  wire  width.  This  happens  because 
it  takes  a  higher  magnetic  field  to  condense  electronic  states  into  cyclotron 
orbits  (Landau  condensation)  when  the  electrostatic  confinement  is  stronger 
(narrower  wires). 

Fig.  8  shows  the  dependence  of  the  normalized  values  of  x^^^  on  wire 
width  for  three  different  values  of  a  magnetic  field  and  a  fixed  value  of  the 
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wave  vector.  For  weak  magnetic  field  of  0.3  tesla,  the  curve  increases 
monotonically  with  increasing  wire  width  (case  “b”  of  Fig.  6).  This  happens 
because  dm,n  is  proportional  to  the  wire  width  W  (see  Equation  5  which  is 
valid  at  zero  field).  The  magnetic  field  is  obviously  not  strong  enough  for 
the  onset  of  Landau  condensation.  At  a  moderate  magnetic  flux  density  of  1 
tesla,  we  can  observe  some  saturation  features  and  for  a  strong  magnetic  field 
of  5  tesla,  the  curve  is  non-monotonic,  rolling  down  to  almost  zero  for  the  wire 
width  of  1000  A  (case  “c”  of  Fig.  6).  The  physics  underlying  the  difference 
in  the  behaviors  of  the  three  curves  is  essentially  the  same  as  that  responsible 
for  the  features  in  Fig.  7.  At  small  values  of  wire  width  {W  «  170A),  there 
is  an  additional  peak  in  the  x^^^  curve.  This  peak  is  a  manifestation  of  the 
fact  that  flat  has  become  comparable  to  the  operational  frequencies.  At  this 
point  the  resonant  behavior  is  strongly  suppressed  by  the  near  vanishing  of 
the  dipole  elements. 

In  our  numerical  calculations  we  have  used  N=  10^^  cm~^.  For  this  dilute 
concentration,  high  density  effects  such  as  screening  and  bandgap  renormal¬ 
ization  are  not  important  and  Equation  (11)  is  quite  valid.  In  fact,  ref.  [1] 
demonstrated  an  excellent  agreement  between  theory  and  experiment  with¬ 
out  accounting  for  any  high  density  effect  even  though  the  carrier  concentra¬ 
tion  in  that  study  was  iV  =  5  x  10^^/cm®.  Therefore,  we  believe  that  high 
density  effects  are  not  significant  in  this  regime. 

The  peak  value  of  the  second  order  susceptibility  for  a  wire  width  of  500 
A  is  x^^^  =  1.5  lO"’^  m/V  (the  absolute  magnitudes  of  the  peak  values  for 
various  wire  widths  are  given  in  the  caption  of  Fig.  7).  For  comparison,  the 
nonlinear  susceptibility  of  electric  field  biased  GaAs  quantum  wells  (W=92 
A)  -  calculated  theoretically  and  measured  experimentally  in  ref.  [1]  -  was 
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=  2.4  10“®  m/V  for  an  electric  field  of  36  kV/cm.  This  shows  that  rela¬ 
tively  weak  magnetic  fields  in  quantum  wires  can  produce  similar  magnitudes 
of  as  rather  strong  electric  fields  in  quantum  wells.  Unfortunately,  there 
is  no  theoretical  or  experimental  result  available  for  either  electric  field  bi¬ 
ased  quantum  wires  or  magnetic  field  biased  quantum  wells  so  that  a  direct 
comparison  is  not  possible.  Nonetheless,  it  is  obvious  that  magnetic  field 
biased  quantum  wires  provide  a  very  attractive  alternative  to  other  means  of 
producing  large  values.  In  fact,  the  largest  value  of  (obtained  at  a 
magnetic  flux  density  of  2  tesla)  in  a  magnetic-field-biased  quantum  wire  is 
three  orders  of  magnitude  higher  than  what  can  be  achieved  in  bulk  GaAs. 

IV.  Conclusion 

We  have  theoretically  studied  the  giant  dipole  effect  in  magnetic-field- 
biased  semiconductor  quantum  wires.  The  dipoles  are  associated  with  tran¬ 
sitions  between  magneto-electric  subbands  within  the  conduction  band,  some 
of  which  are  forbidden  in  the  absence  of  the  magnetic  field.  The  resonant 
frequencies  of  these  transitions  can  be  tuned  by  the  magnetic  field  which 
allows  the  realization  of  externally  tunable  inter-subband  lasers.  We  have 
also  studied  the  possibility  of  second  harmonic  generation  in  a  quantum  wire 
biased  with  a  magnetic  field  and  find  a  strong  second  harmonic  component 
of  the  susceptibility  which  has  applications  in  non-linear  optics. 
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FIGURE  CAPTIONS 


[1]  An  electron  waveguide  (quantum  wire)  subjected  to  a  magnetic  field  along 
the  ^  axis.  The  width  of  the  wire  is  much  larger  than  the  thicknes. 

[2]  Dipole  moments  for  various  inter-subband  transitions  as  functions  of  the 
propagating  wave  vector  k  for  a  magnetic  flux  density  of  1  tesla.  At  zero 
translational  velocity  (A;  =  0)  the  dipole  of  the  transition  eTe3  vanishes. 
The  GaAs  quantum  wire  is  1000  A  wide. 

[3]  The  dipol  es  of  three  inter-subband  transitions  as  functions  of  the  applied 
magnetic  field.  The  dipole  dei-es  peaks  at  a  magnetic  flux  density  of  0.3 
tesla.  The  wire  width  is  the  same  as  that  in  Fig.  2. 

[4]  Energy  vs.  wave  vector  relation  of  electrons  in  (a)  first  subband,  and  (b) 
third  subband  of  a  1000  A  wide  quantum  wire.  The  wave  vector  is  along 
the  free  propagation  direction.  The  results  are  shown  for  three  values  of  a 
magnetic  field.  The  energy  is  calculated  from  the  bulk  conduction  band  edge. 

[5]  The  y  component  of  the  electron  envelope  function  for  the  first  subband 
at  a  magnetic  flux  density  of  0  and  10  tesla. 

[6]  The  y  component  of  the  electron  envelope  functions  for  the  first  and 
third  electronic  subbands.  The  results  are  shown  for  the  cases  when  (a)  no 
magnetic  field  is  present,  (b)  when  a  weak  magnetic  field  is  present  and, 
finally,  (c)  when  a  strong  magnetic  field  is  present. 

[7]  Second  order  susceptibility  as  a  function  of  the  biasing  magnetic  field. 
The  peak  values  of  the  susceptibility  are  13.2  10~^  m/V,  1.5  10“^  m/V  and 
3  10“®  m/V  for  wire  widths  of  1000  A,  500  A  and  300  A  respectively.  The 
results  are  shown  for  the  wave  vector  A;=0.01  /  A  (fixed  excitation  frequency). 
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[8]  Second  order  susceptibility  as  a  function  of  the  wire  width  for  three  values 
of  the  biasing  magnetic  field.  The  maximum  values  of  the  curves  are  the 
same  as  in  Fig.  7.  The  narrow  peaks  at  a  wire  width  of  200  A  are  due  to 
resonances  occuring  when  either  fija  =  +a;2  or  flea  =  ^2  (see  Equation  11 

in  the  text). 
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Abstract 

In  this  paper,  we  have  calculated  the  refractive  index  of  a  quantum 
wire  waveguide  in  the  vicinity  of  an  exciton-polariton  resonance.  The 
critical  values  of  the  exciton  decay  parameter,  defining  the  onset  of 
polariton  transport  regime,  and  the  associated  temperature  were  also 
found  Tising  the  combination  of  a  variational  approach  and  a  numeri¬ 
cal  solution.  Our  theoretical  model  allows  us  to  include  the  effects  of 
an  external  magnetic  field.  The  results  show  that  confinement  of  exci- 
tons  to  one  dimension  and  the  simultaneous  application  of  a  magnetic 
field  may  lead  to  the  extension  of  the  temperature  and  spatial  limits 
of  polariton  transport.  The  magnetic  field  can  be  used  to  shift  refrac¬ 
tive  index  peaks  in  frequency  thus  providing  a  much-desired  tuning 
capability. 


I.  Introduction 

Exciton  dynamics  in  semiconductor  quantum  confined  structures  has  al¬ 
ways  been  attractive  because  of  its  potential  applications  in  optoelectronics. 
One  of  the  most  interesting  phenomenon  related  to  excitons  in  such  struc¬ 
tures  is  the  formation  of  exciton-polaritons.  In  the  spectral  region  around  an 
exciton  resonance,  a  photon,  absorbed  by  a  semiconductor,  linearly  couples 
with  an  exciton  to  create  a  polariton.  Exciton  polaritons  have  been  studied 
extensively  in  quantum  wells  by  measuring  photoluminescence  and  reflection, 
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as  well  as  by  picosecond  time-of-flight  measurements.^’^  It  was  shown  that 
polaritons  are  much  more  stable  in  quantum  wells  than  in  bulk,^  and  they  are 
expected  to  be  even  more  stable  in  quantum  wires  since  the  exciton  binding 
energy  and  oscillator  strength  tend  to  increase  with  reducing  dimensionality.^ 

The  formation  of  polaritons  modifies  the  transport  of  light  through  the 
medium.  In  particular,  the  medium  of  propagation  becomes  substantially 
more  transparent®’®  and  the  group  velocity  of  light  propagating  along  the 
waveguide  approaches  the  speed  of  light  in  vacuum.^  At  the  same  time,  there 
is  a  possibility  of  controlling  polariton  transport  with  an  external  field.  In  this 
paper,  we  argue  that  a  relatively  weak  magnetic  field  is  particularly  attractive 
for  tuning  polariton  transport  because  it  increases  exciton  oscillator  strength 
thus  extending  the  polariton  regime  of  energy  transfer,  while  an  electric  field 
would  decrease  exciton  oscillator  strength  and  quench  polariton  transport. 

Although  existing  theoretical  models  recognize  modification  of  polariton 
transport  due  to  spatial  confinement,  they  account  for  it  by  using  heuristi- 
cally  peaked  values  for  exciton  oscillator  strength  and  binding  energy.  The 
authors  are  not  aware  of  any  attempts  to  include  an  external  field,  particu¬ 
larly  magnetic,  into  consideration. 

In  this  work,  we  present  a  model  for  calculating  the  exciton-polariton  crit¬ 
ical  decay  parameter  and  the  refractive  index  of  a  quantum  wire  around  a 
polariton  resonance  in  the  presence  of  a  magnetic  field.  The  decay  parameter 
determines  the  regime  of  polariton  transport.  To  our  knowledge,  this  is  the 
first  study  where  the  exciton  longitudinal-transverse  (LT)  splitting  and  ex¬ 
citon  resonance  frequency  -  which  define  the  polariton  dispersion  -  are  found 
in  a  non-ad-hoc  manner  using  the  combination  of  a  variational  approach  and 
an  exact  numerical  solution  of  the  Schrodinger  equation.  The  calculations 
are  performed  for  a  quantum  wire  with  finite  lateral  dimensions  subjected  to 
a  magnetic  field. 

The  rest  of  the  paper  is  organized  as  follows.  In  section  II,  we  establish  the 
polariton  dispersion  relation  used  throughout  the  model;  section  III  presents 
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the  results  of  calculation  of  the  LT  splitting  and  oscillator  strength  of  the 
exciton  transition  in  a  quantum  wire  subjected  to  a  magnetic  field;  in  section 
IV,  we  examine  the  exciton  critical  damping  for  the  onset  of  the  polariton 
transport  regime  and  calculate  refractive  index  of  the  wire  in  the  vicinity  of 
exciton  resonance.  Conclusions  are  given  in  section  V  of  the  paper. 

II.  Polariton  dispersion 

We  consider  an  array  of  parallel  GaAs  quantum  wires  of  rectangular  cross 
section  separated  by  infinite  potential  barriers  so  that  wavefunctions  of  the 
excitons  from  different  wires  do  not  overlap.  In  such  a  structure,  excitons 
are  free  to  move  along  the  wire  axes  but  are  confined  in  perpendicular  direc¬ 
tions.  The  lateral  dimension  of  each  wire  is  comparable  to  the  exciton  Bohr 
radius.  The  cladding  material  is  assumed  to  have  a  similar  refractive  index 
so  that  we  can  ignore  image  charges  of  the  exciton  and  associated  dielectric 
confinement  effects.  Under  these  conditions,  the  dispersion  relation  of  the 
exciton  polaritons  can  be  determined  for  each  separate  wire.  The  multiple 
wire  structure  in  this  case  merely  forms  a  waveguide  structure  analogous  to 
that  considered  in  Ref.  [8]. 

Most  theoretical  models  for  exciton  polaritons  (both  in  bulk  material 
and  nanostructures)  embody  a  semi-classical  approach  and  utilize  the  disper¬ 
sion  relation  of  a  polariton  derived  for  a  single  electric-dipole-active  exciton 
resonance.®’®  Here,  we  adopt  the  same  philosophy  and  consider  electromag¬ 
netic  waves  propagating  through  an  array  of  quantum  wires  with  a  wavevec- 
tor  k  parallel  to  the  wire  axis.  This  choice  of  the  direction  of  propagation 
allows  for  a  spatial  dispersion  of  the  light  waves.  In  the  opposite  case  of  light 
propagating  normal  to  the  wire  axis,  the  translational  motion  of  excitons  is 
suppressed  and  the  spatial  dispersion  effects  do  not  occur. 

In  the  long- wave  approximation  {kLy^z  <  1,  where  Ly^z  are  the  wire  lateral 
dimensions)  the  array  interacts  with  light  waves  like  an  effective  medium,  and 
the  dielectric  function  in  the  vicinity  of  an  isolated  exciton  resonance  can  be 
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written  as^ 

/  l\  I  2eoOJLT^o  ^ 

e{u),k)—eo-\ — 5 - 2  1  *72 — TUI — ^ — •  (1) 

o;“  —  +  hk^Uo/M  —  zuj 

where  Co  is  the  background  dielectric  constant  (contribution  made  by  other 
resonances),  cu  is  the  frequency  of  light,  ult  is  the  longitudinal-transverse 
splitting  of  the  exciton  related  to  its  oscillator  strength,  Uo  is  the  exciton 
resonant  frequency  &t  k  =  0,  M  =  mg  +  m/i  is  the  translational  mass 
of  an  exciton,  and  T  =  hj  is  the  exciton  damping  parameter.  Here  we 
have  assumed  parabolic  wavevector  dependence  of  the  exciton  frequency 
fujjt{k)  =  tkjJo  +  h‘^k^/2M,  with  the  caveat  that  this  is  valid  only  in  weak 
magnetic  fields  when  the  magnetic  length  Im  (=  ^h/eB)  is  much  larger  than 
the  transverse  dimensions  of  the  wire.  In  the  formula  above,  the  quanti¬ 
ties  (jJlt  =  <^LT{Ly,z,  B)  and  Uo  =  oJo{Ly,z,B)  are  the  functions  of  the  wire 
lateral  dimensions  and  a  magnetic  field.  The  exciton  damping  constant  is 
considered  to  be  independent  of  the  magnetic  field  since  it  is  known  that 
energy- averaged  phonon-interaction  rates  in  quantum  wires  are  not  terribly 
sensitive  to  a  magnetic  field.  In  any  case,  Eq.  (1)  is  a  good  approximation 
when  the  magnetic  field  applied  to  the  system  is  relatively  weak:  >  Ly^z- 

This  equation  relates  uj  and  k  and  is  the  sought-after  dispersion  relation  of 
a  polariton. 

Before  we  can  go  further  into  polariton  transport  properties,  we  have  to 
calculate  ojit  and  Uo  as  the  functions  of  wire  dimensions  Ly^^  and  a  magnetic 
flux  density  B.  This  is  discussed  in  the  next  section. 


III.  Longitudinal-transverse  splitting 


Let  us  assume  that  the  infinite  potential  barriers  of  the  quantum  wire  are 
located  at  ^  =  dcLy/2  and  z  =  d=Z-^/2.  A  magnetic  field  is  applied  along 
the  z-direction  (see  inset  to  Fig.  1).  To  simplify  the  calculations,  we  assume 
strong  quantum  confinement  of  the  carriers  which  enables  us  to  factorize  an 
exciton  wavefunction  into  the  product  of  electron  and  hole  wave  functions. 
Moreover,  we  limit  our  consideration  to  systems  with  relatively  large  dielec- 
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trie  constants  so  that  all  Coulomb  interactions  are  strongly  screened.  This 
assumption,  together  with  the  hard-wall  boundary  condition,  allows  a  co¬ 
ordinate  separation.  Consequently,  the  wave  function  of  an  exciton  in  the 
vicinity  of  subband  bottom  (with  center-of-mass  momentum  Px  ~  0)  is  given 
by^ 

^  Ve,  Vh,  Ze,  Zh)  =  gt{x,  v)A{ye,  Zeji^hiVh,  Zh)  (2) 

=  9t{x,  ri)(l>e{ye)(f>h{yh)Xe{Ze)Xh{Zh), 
where  gt{x,r])  is  chosen  to  be  the  Gaussian- type  “orbital”  function: 

=  (3) 

TT 

in  which  77  is  a  variational  parameter  which  defines  the  exciton  size  (“longitu¬ 
dinal  length”),  and  x  is  the  relative  electron-hole  coordinate.  The  subscripts 
in  Xe,h,ye,h,Ze,h  identify  them  as  electron  or  hole  coordinates.  The  variables 
Xe,h{ze,h)  are  the  z-components  of  the  wave  functions  which  are  not  affected 
by  the  magnetic  field.  They  are  given  by  particle-in-a-box  states.  The  elec¬ 
tron  and  hole  wave  functions  along  the  y  direction,  4>e^{ye,h),  are  to  be  cal¬ 
culated  numerically  when  a  magnetic  field  is  present.  This  is  done  by  solving 
the  one-particle  Schrodinger  equation  using  a  finite  difference  method. 

In  order  to  find  an  exciton  “length”  77,  we  use  the  variational  approach  of 
minimizing  the  energy  given  by  <  >,  where  the  exciton  Hamiltonian 

is 

ux  —  4.  4.  Py‘  4.  Pyi'  ('4'j 

2M  2n  2m^  2mh  ^  ^ 

+  +  eB(ye/me  -I-  yh/mh)Px  +  ^^(yeMe  +  yV'^h) 

+  Uc{Xe,  Xh,  ye,  yh,  Ze,  Zh)  +  Us{ye,  yh,  Ze,  Zh)- 

Here  we  have  chosen  the  Landau  gauge  A  =  (—By,  0,0).  The  quanti¬ 
ties  me.,mh,  are  the  effective  masses  of  electrons  and  holes  respectively, 
!///(=  l/rue  +  1/m/i)  is  the  exciton’s  reduced  mass,  Uc{xe,Xh,ye,yh,Ze,Zh) 
is  the  electron-hole  Coulomb  interaction  term,  Us{ye,yh,Ze,Zh)  is  the  spa¬ 
tial  confinement  potentials  for  electrons  and  holes  along  y  and  z  directions. 
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Details  of  the  variational  procedure,  calculations  of  77,  etc.,  can  be  found  in 
some  of  our  earlier  work."^ 

We  can  find  the  oscillator  strength  of  the  exciton  transition  ag  and  the 
LT  splitting  by  evaluating  the  momentum  matrix  element  which  is  given  as 

=  1^  J  dkgt{x,r))Mc„{k)\‘^,  (5) 

where  Mc„  is  the  valence-band  to  conduction-band  dipole  matrix  element, 
and  k  is  again  the  wave  vector  along  the  unconfined  direction  of  the  wire. 
When  the  k  dependence  of  Ma,  is  neglected,  Eq.  (5)  reduces  to  the  simple 
expression 

=  =  (6) 


The  exciton  oscillator  strength  per  unit  length  can  be  written  as  follows 


motUOn 


Here  huo  =  Eq  +  Ee\  -f  E^hi  —  min<  >  is  the  exciton  ground  state 

energy,  Eq  is  the  fundamental  bandgap  of  the  bulk  material,  Egi,  Ehhi  are 
the  lowest  electron  and  the  highest  heavy  hole  magneto-electric  subband 
bottom  energies  in  a  quantum  wire  measured  from  the  bottom  of  the  bulk 
conduction  band  and  the  top  of  the  bulk  valence  band,  and  is  the  free 
electron  mass.  The  exciton  LT  splitting  Clir  =  ^lt  can  now  be  written  as 


In  Fig.  1  we  present  the  LT  splitting  calculated  for  different  wire  dimen¬ 
sions  and  magnetic  flux  densities.  The  physical  parameters  used  for  the  calcu¬ 
lations  correspond  to  a  GaAs  quantum  wire  with  e  =  12.9eo,  Eg  —  1.515eV, 
rUe  =  O.OGTrrZo,  rrih  =  0.5mo,  where  rUo  is  free  electron  mass  and  Co  is  electri¬ 
cal  permitivitty  of  free  space,  Egi,  Ehhi  are  calculated  numerically  following 
the  prescriptions  of  Refs.  [10,  11].  One  can  see  from  the  figure,  that  the  ex¬ 
citon  splitting  is  sensitive  to  the  spatial  confinement  and  increases  by  about 
60  %  when  the  wire  width  decreases  from  500  A  to  50  A.  A  magnetic  field 
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Figure  1:  Longitudinal-transverse  splitting  of  an  exciton  as  a  function  of 
wire  width  (left  pannel);  and  as  a  function  of  magnetic  flux  density  for  a 
500  A  wide  wire  (right  panel).  For  both  panels,  the  upper  and  lower  curves 
correspond  to  200  A  and  500A  thick  wires,  respectively. 

also  increases  the  splitting  (and  oscillator  strength)  thus  making  the  exciton 
polaritons  more  stable. 

IV.  Refractive  index  and  decay  parameter 

The  experimentally  observed  higher  transparency  of  the  medium  of  prop¬ 
agation  in  polariton  transport  regime  has  been  attributed  to  certain  features 
of  the  dispersion  law  for  excitonic  polaritons  and  to  the  fact  that  polariton 
transport  by  itself  cannot  cause  true  absorption.  In  order  for  absorption  to 
occur,  polaritons  have  to  be  scattered  inelastically,  e.g.,  by  phonons.  The  on¬ 
set  of  polariton  transport  through  some  structure  is  governed  by  the  exciton 
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polariton  coherence  length  related  to  the  exciton  decay  parameter  F  =  h'y. 
It  has  been  shown,  both  experimentally  and  theoretically,  that  there  exists 
a  critical  value  of  the  exciton  decay  parameter,  Fc,  which  corresponds  to  a 
change  in  the  nature  of  absorption.®’®  Here  we  intend  to  examine  the  influence 
of  spatial  confinement  and  a  magnetic  field  on  this  parameter  and  calculate 
the  refractive  index  of  a  quantum  wire  in  the  polariton  regime. 

Confining  ourselves  to  the  TE  waves  relevant  to  light  propagation  in  the 
medium,  we  can  write  polariton  dispersion  in  the  following  form 

eioj,  k)  =  ^  =  n^.  (9) 

(jJ 


Combining  this  equation  with  Eq.  (1),  and  after  some  algebra,  we  obtain 

+  (^2 _^2 _ -i'^U}^-2u}LT^o)  =0.  (10) 


This  equation  can  be  solved  for  two  sets  of  the  refractive  indices,  ni  and  na, 
corresponding  to  different  transverse  polariton  branches.  It  also  follows  from 
Eq.  (10)  that  if  the  damping  parameter  F  becomes  larger  than  the  critical 
value 

r„  =  B)  =  (11) 

then  only  one  light  wave  mode  can  propagate  in  the  medium,  since  there  is 
only  one  real  solution  for  n.  This  is  the  boundary  of  the  polariton  propaga¬ 
tion  regime.  The  critical  value  comes  about  because  of  the  term  hk^Uo/M 
associated  with  spatial  dispersion  effects.  The  physical  importance  of  the 
critical  damping  can  be  illustrated  by  the  following  example.  It  was  shown® 
that  when  the  damping  exceeds  the  critical  value,  the  integral  absorption  is 
independent  of  F  and  proportional  to  the  oscillator  strength  of  transitions 
(non-polariton  regime) .  When  F  <  Fc,  the  integral  absorption  depends  on  F 
linearly,  and  decreases  with  decreasing  damping. 

Using  the  results  from  the  previous  section,  we  can  calculate  Fc  for  dif¬ 
ferent  values  of  wire  widths  and  magnetic  field.  In  Fig.  2  we  present  the 
critical  exciton  damping  (decay)  parameter  as  a  function  of  wire  width.  It  is 
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Figure  2:  The  critical  exciton  decay  parameter  (damping)  as  a  function  of 
wire  width.  The  upper  and  lower  curves  correspond  to  the  thickness  along 
the  z-direction  of  200  A  and  300  A,  respectively. 

normalized  by  the  value  of  the  decay  parameter  of  a  very  wide  wire  (Lj,=700 
A)  which  is  wide  enough  to  be  approaching  the  2D  limit.  For  this  wide 
wire,  Fc  «  1.01  meV.  Combining  Eqs.  (7,  8,  11)  we  can  also  estimate  the 
magnetic  field  dependence  of  the  critical  decay  parameter  using  the  formula 
rc(jB)/rc(0)  =  Although  not  shown  here, 

the  magnetic  field  dependence  of  the  critical  parameter  is  weak;  it  increases 
only  5%  at  a  magnetic  flux  density  of  5  tesla.  The  strong  dependence  of  the 
critical  damping  on  the  wire  width  may  lead  to  a  pronounced  modification  of 
the  integral  (total)  absorption  of  systems  consisting  of  a  number  of  narrow 
quantum  wires. 

We  can  now  find  a  temperature  which  corresponds  to  the  critical  damp¬ 
ing  from  the  relation  Fg  =  Fj,  -f-  Fp/i(Tc),  where  Fq  is  the  damping  associ¬ 
ated  with  the  impurity  and  other  temperature-independent  elastic  scatter¬ 
ing,  while  TphlTc)  represents  interactions  with  acoustic  and  optical  phonons. 
By  increasing  Fc  one  can  increase  Tc  which  defines  the  onset  of  polariton 
transport  and,  as  a  consequence,  higher  transparency. 
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Since  there  is  no  data  available  on  quantum  wires,  we  assume  that  the 
half- width  at  half  maximum  (HWHM)  of  the  exciton  resonance  in  a  quantum 
wire  is  the  same  as  in  a  200  A  thick  GaAs/AlGaAs  quantum  well.  Using  the 
approximation  of  Ref  [11]  we  may  write  for  our  case  (energy  units  are  meV) 

Te  =  r+  -f  0.00147Tc  -1-  4.0(6''"^'*/*=^^'  -  1)"^  (12) 

where  huph  =  36  meV  is  a  longitudinal  optical  phonon  energy,  ks  is  the 
Boltzman  constant,  Eb  «  10  meV  is  the  everage  binding  energy  for  donor 
impurities  in  GaAs,  Timp  =  0.75  meV  is  a  linewidth  due  to  fully  ionized 
impurity  scattering,  r+  «  0.45  meV  is  the  linewidth  due  to  inhomogeneous 
fluctuations  of  the  wire  thickness.  The  values  chosen  for  the  various  param¬ 
eters  are  typical  of  experimental  systems  reported  in  the  literature. 

Table  I.  Critical  temperature  vs.  wire  width 


L„  (A) 

500 

300 

100 

Fc,  (meV) 

1.09 

1.16 

1.45 

Z,  (K) 

118 

138 

169 

Solving  Eq.  (12)  for  the  temperature  Tc  for  each  value  of  TdLy^z,  B),  we 
are  able  to  obtain  the  dependence  of  the  critical  temperature  on  the  wire 
width.  The  thickness  of  the  wire,  was  fixed  at  200A  for  this  calculation. 
As  one  can  see  from  Table  I,  that  the  critical  temperature  Tc,  that  defines  the 
onset  of  polariton  transport,  can  be  controlled  over  a  wide  range  by  changing 
the  wire  width  Ly. 

Now  let  us  assume  that  F  <  Fc  (exciton  polariton  regime)  and  find  the 
refractive  index  of  the  quantum  wire  in  the  vicinity  of  polariton  resonance. 
In  order  to  do  this,  we  make  use  of  Pekar’s  additional  boundary  condition 
(total  polarization  is  zero  at  the  boundary)  and  write  the  effective  refractive 
index  as 


Tleff  — 


niTia  +  Cq 

ni  +  n2 


(13) 


In  Fig.  3,  we  present  the  real  (upper  panel)  and  imaginary  (lower  panel) 
parts  of  the  refractive  index  of  a  quantum  wire  with  thickness  200  A  and 
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Energy  (meV) 


Figure  3:  The  refractive  index  of  the  wire  in  the  vicinity  of  exciton  resonance. 
The  solid,  dashed  and  dash-dotted  curves  correspond  to  a  0,  1,  and  5  tesla 
magnetic  flux  density,  respectively. 

width  300  A.  An  external  transverse  magnetic  field  is  applied  along  the 
thickness.  The  decay  parameter  is  chosen  to  be  F  =  ImeV  which  is  less 
than  Fc  for  the  given  wire  dimensions.  As  one  can  see  from  the  figure,  the 
maximum  of  the  real  part  of  the  refractive  index  is  as  large  as  4.97  at  zero 
field  and  5.04  at  5  tesla  magnetic  field.  It  is  about  1.4  times  larger  than 
that  of  the  bulk  material.  The  refractive  index  attains  its  maximum  value 
at  a  photon  energy  slightly  lower  than  the  resonance  energy  fkjJo{Ly^z,B) 
at  any  given  magnetic  field  B.  The  minimum  value  of  the  refractive  index, 
which  is  1.66  for  zero  field  and  1.59  for  5  tesla,  is  located  at  a  frequency  of 
^o{Ly,z,  B)  •\-(jJLT{.Ly^z,  B).  It  is  interesting  to  note  that  efficient  waveguiding 
can  be  achieved  in  the  spectral  range  where  the  real  part  of  the  refractive 
index  increases.  However,  in  this  region,  the  imaginary  part  of  the  refractive 
index  (extinction  coefficient)  also  peaks  and  this  increases  the  transmission 
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loss.  The  mitigating  factor  in  all  this  is  that  the  two  peaks  do  not  occur  at 
exactly  the  same  frequency  so  that  an  optimal  region  for  optical  waveguiding 
exists. 

A  magnetic  field  of  5  tesla  blue-shifts  the  refractive  index  peak  by  2  meV. 
According  to  Ref.  [13],  approximately  the  same  magnitude  of  an  opposite 
red-shift  can  be  achieved  by  applying  an  electric  field  of  about  4x10^  V/cm. 
However,  the  electric  field  leads  to  a  15-20  %  increase  in  exciton  radius  and  a 
concomitant  decrease  in  the  binding  energy.  This,  in  turn,  causes  a  decrease 
in  exciton  LT  splitting  and  makes  exciton  polaritons  less  stable.  Note  that 
while  an  electric  field  will  tend  to  ionize  an  exciton  by  pulling  the  electron 
and  hole  apart,  a  magnetic  field  has  the  opposite  effect.  It  squeezes  the 
electron  and  hole  even  tighter  together  and  increases  the  binding  energy. 
Therefore,  the  magnetic  field  can  be  used  to  advantage  in  this  context  since 
it  shifts  the  peaks  in  frequency  while  actually  increasing  polariton  stability. 
This  frequency  tuning  capability,  acquired  without  a  penalty  in  polariton 
stability,  is  obviously  very  attractive  and  has  device  applications. 

VI.  Conclusions 

In  this  paper,  we  have  calculated  the  refractive  index  of  a  quantum  wire 
waveguide  in  the  vicinity  of  polariton  resonance.  The  critical  values  of  the 
exciton  decay  parameter  and  associated  temperature  were  also  found  taking 
into  account  the  effects  of  spatial  confinement  and  an  external  magnetic 
field.  Our  results  show  that  confinement  of  excitons  to  one  dimension  and 
the  application  of  a  magnetic  field  may  lead  to  the  extension  of  temperature 
and  spatial  limits  of  polariton  transport.  The  magnetic  field  can  be  used  to 
shift  refractive  index  peaks  in  frequency  -  without  compromising  polariton 
stability  -  thus  providing  a  much-desired  tuning  capability. 
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Abstract 

We  have  theoretically  studied  non-lineair  frequency  conversion  in 
a  semiconductor  quantum  wire  biased  with  a  magnetic  field.  In  these 
systems,  eflScient  second  harmonic  generation  occurs  as  a  result  of 
the  large  value  of  the  second-order  dielectric  susceptibility  aris¬ 
ing  from  dipole  transitions  between  magneto-electric  subbands.  The 
magnitude  and  peak  frequency  of  as  well  as  the  absorption  and 
refractive  index  associated  with  can  be  tuned  with  the  magnetic 
field.  This  allows  one  to  achieve  low  insertion  loss  and  efficient  phase 
matching  by  manipulating  the  absorption  and  refractive  index  with  a 
magnetic  field. 


I.  Introduction 

Most  ordinary  solids  are  not  efficient  frequency  converters  because  they 
exhibit  extremely  small  even-order  dielectric  susceptibilties.  Ideally,  even- 
order  susceptibilities  vanish  in  solids  with  inversion  symmetry.^  Conse¬ 
quently,  a  semiconductor  structure  can  exhibit  a  large  value  of  the  second 
order  susceptibility  x^^^  only  if  the  inversion  symmetry  of  the  conduction- 
band  potential  is  broken  artificially  either  by  an  external  electric  field,  or  by 
the  intentional  growth  of  an  asymmetric  structure.  Obviously,  the  former  is 
the  preferred  method  since  an  electric  field  can  be  varied  continuously  and 
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this  allows  one  to  tune  the  degree  of  symmetry-breaking  and  the  magnitude 
of  However,  an  electric  field  has  a  practical  shortcoming.  In  a  quantum 
confined  structure,  it  tilts  the  potential  barriers  confining  the  photogener¬ 
ated  carriers.  As  a  result,  carriers  can  escape  by  tunneling  or  thermionic 
emission  and  this  is  especially  serious  in  GaAs/AlGaAs  systems  where  the 
barrier  height  is  relatively  small.  Indeed,  it  has  been  pointed  out  that  the 
electronic  states  in  a  quantum  confined  system  biased  by  a  transverse  elec¬ 
tric  field  are  never  true  bound  states  since  the  particles  can  always  lower 
their  energy  by  escaping  from  the  well^.  Consequently,  optical  transitions 
(and  their  higher  order  harmonics)  associated  with  these  states  have  incon¬ 
veniently  large  linewidths  and  small  oscillator  strengths. 

Recently,  we  proposed  magnetostatic  biasing  as  an  attractive  alternative 
to  mitigate  this  problem.^  We  showed  that  a  magnetic  field  can  break  in¬ 
version  symmetry  in  a  quantum  wire  without  tilting  potential  barriers.  A 
transverse  magnetic  field,  applied  to  a  wire,  exerts  a  Lorentz  force  on  an 
electron  moving  along  the  length.  As  a  result,  its  wave  function  (in  any 
magneto-electric  subband)  will  be  skewed  towards  one  edge  of  the  wire.  This 
skewing  does  not  tilt  potential  barriers  to  first  order  (the  barriers  may  tilt 
slightly  because  of  a  second-order  effect  associated  with  space-charges  and 
the  self-consistent  (Hall)  electric  field).  However,  it  effectively  breaks  in¬ 
version  symmetry  since  it  causes  net  charges  to  accumulate  at  either  edge 
of  the  wire.  This  leads  to  a  non-vanishing  even-order  susceptibility  in  an 
otherwise  symmetric  structure.  The  skewing  has  another  subtle  effect.  The 
degree  to  which  the  wave  function  is  skewed  is  different  in  different  subbands 
since  an  electron  has  different  kinetic  energies  (and  hence  experiences  differ¬ 
ent  Lorentz  forces)  in  different  subbands.  As  a  result,  transitions  between 
subbands  whose  wave  functions  have  the  same  parity  -  which  are  forbidden 
without  a  magnetic  field  -  are  now  allowed  since  the  parities  are  altered  by 
different  amounts  in  different  subbands.  This  effect  has  some  similarity  with 
the  quantum  confined  Lorentz  effect  (QCLE)  previously  examined  by  us^  in 
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the  context  of  interband  transitions  between  conduction  and  valence  band 
states. 

In  this  paper,  we  first  calculate  the  second-order  susceptibility  in  a 
symmetric  quantum  wire  whose  inversion  symmetry  (along  the  width)  has 
been  broken  with  a  magnetic  field.  We  restrict  ourselves  to  narrow  GaAs 
wires  with  a  width  of  about  150 A.  The  energy  spacing  between  the  first  and 
second  subband  is  AE12  ~  72  meV.  This  choice  of  the  wire  dimension  puts 
the  resonant  frequency  (for  transitions  between  the  lowest  subbands)  in  the 
mid-infrared  spectral  region.  The  wavelength  of  the  second  harmonic  compo¬ 
nent  of  this  transition  is  about  8.6  ^m.  Here  we  will  be  mainly  interested  in 
arising  from  resonant  and  near-resonant  inter-subband  transitions  which 
are  governed  by  the  interplay  of  dipoles  and  resonant  excitations.  In  con¬ 
trast,  Ref.  [3]  focussed  on  the  off-resonance  regime  which  was  governed  solely 
by  the  dipoles.  We  will  also  calculate  absorption  and  refractive  index  in  the 
frequency  region  of  interest  for  both  pump  and  second  harmonic  frequencies 
and  show  how  they  can  be  manipulated  with  an  external  magnetic  field  to 
realize  low  insertion  loss  and  efficient  phase  matching. 

The  rest  of  the  paper  is  organized  as  follows.  In  the  next  section,  we 
describe  the  theoretical  formulation,  followed  by  results.  Finally,  in  section 
IV,  we  present  the  conclusions. 

II.  Theory 

We  consider  a  generic  GaAs  quantum  wire  (as  shown  in  the  inset  of  Fig. 
1)  with  a  magnetic  field  applied  along  the  z  direction.  The  thickness  along 
the  z  direction  is  so  small  (and  consequently  the  subband  separation  in  energy 
in  this  direction  is  so  large)  that  for  the  range  of  photon  energies  considered, 
an  electron  cannot  be  excited  (by  real  transition)  into  a  subband  which  has 
more  than  two  nodes  along  the  z-direction.  In  other  words,  such  a  transition 
will  not  be  accessible  in  energy.  This  restriction,  coupled  with  the  fact  that  a 
magnetic  field  does  not  affect  the  z-component  of  the  electron  wave  function, 
allows  us  to  drop  the  z-component  from  further  consideration.  The  width 
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of  the  wire  along  the  y-direction  is  however  large  enough  (W  =  150 A)  that 
subbands  with  more  than  two  nodes  along  the  y-direction  are  accessible  in 
energy. 

In  systems  without  inversion  symmetry,  the  lowest  order  optical  nonlin¬ 
earity  is  of  the  second  order  and  is  expressed  by 

where  P  is  the  polarization  caused  by  two  electric  fields  Ei  and  E2  that 
are  associated  with  the  electromagnetic  fields  of  either  two  frequency  compo¬ 
nents  of  the  same  light  beam  or  two  different  coherent  beams  with  frequencies 
Ui  and  wave  vectors  ki.  It  is  well  known  that  the  third-ranked  tensor 
will  vanish  in  any  structure  with  inversion  symmetry.  A  quantum  confined 
structure  may  lack  inversion  symmetry  for  two  main  reasons,  (i)  the  semh 
conductor  material  by  its  intrinsic  chemical  and  crystalline  structure  may 
lack  inversion  symmetry,®  and  this  is  the  case  in  most  III-V,  II- VI,  and  I- 
VII  compounds  along  certain  crystallographic  directions,  or  (ii)  the  quantum 
confining  potential  well  may  be  asymmetric  (e.  g.  triangular  potential  well, 
asymmetric  double  square  well  potential,  etc.).  In  the  first  case,  the  asym¬ 
metry  is  related  to  the  intracell  charge  asymmetry  and  is  not  affected  by  the 
confinement  since  the  latter  extends  over  several  unit  cells.  In  the  second 
case,  the  asymmetry  is  artificially  imposed  and  therefore  can  be  engineered. 
It  clearly  depends  on  the  confining  potential.  Insofar  as  an  applied  electric 
field  can  alter  the  potential,  it  can  change  the  degree  of  symmetry-breaking 
and  hence  modulate  x^^^- 

In  the  present  work  we  restrict  ourselves  to  the  second  case  and  do  not 
consider  intrinsic  second  order  nonlinearities  which  can  be  quite  large  in  some 
materials  (the  nonlinear  susceptibility  of  bulk  GaAs  is  xi? =3.8  m/V).® 

As  mentioned  before,  we  avoid  the  use  of  a  symmetry-breaking  electric  field 
since  it  promotes  carrier  escape.  Instead,  we  consider  a  magnetic  field.  Al¬ 
though  a  magnetic  field  does  not  directly  affect  the  potential,  it  leads  to  an 
uneven  charge  distribution  along  the  width  {y  -  axis)  of  the  wire  because  of 
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the  different  degrees  of  skewing  of  the  wave  functions  in  different  magneto¬ 
electric  subbands.  This  has  the  effect  of  breaking  inversion  symmetry. 

As  mentioned  before,  a  magnetic  field  induces  forbidden  transitions  be¬ 
tween  subbands  of  the  same  parity.  The  large  magnitude  of  the  dipoles 
associated  with  these  transitions  and  their  extreme  sensitivity  to  the  field 
open  up  the  possibility  of  controllable  second  harmonic  generation  (SHG) 
that  can  be  manipulated  by  the  magnetic  field.  In  order  to  evaluate  the 
magnitude  and  dependence  of  SHG  on  the  biasing  field  and  wire  geometry, 
we  calculate  the  second  order  susceptibility  using  the  formula^ 


X^^\2uj]u}-,u}) 


a  'ST' _ 


ba-^-  ii)  (f^ca  -  2a;  -  27)  ’ 


(2) 


where  N  is  concentration  (number  density)  of  conduction  electrons,  7  is 
a  damping  factor  associated  with  elastic  and  inelastic  scattering,  = 
hQtja{B,  W,  k)  is  the  energy  spacing  between  the  b-th  and  a-th  magnetoelec¬ 
tric  subbands  which  depends  on  the  applied  magnetic  field,  wire  width  and 
electron  wave  vector  and  dmn  =  dmn{B,W,k)  is  a  dipole  element  of  tran¬ 
sitions  between  different  subbands.  The  total  symmetrisation  operation  St 
indicates  that  the  expression  which  follows  is  to  be  summed  over  all  permuta¬ 
tions  of  the  pairs  (fi,  2u),  (cx,co),  (j3,u>).  Since  St  involves  a  summation  over 
all  possible  permutations,  it  is  clear  that  xjf^^(2a;;a;;a;)  is  invariant  under 
any  of  them. 

In  order  to  calculate  dipole  elements  dmn{B,  W) ,  we  proceed  as  in  Ref.  [3]. 
Under  the  electric  dipole  approximation,  the  matrix  element  of  photoinduced 
inter-subband  transitions  within  the  conduction  band  is  given  by^ 


df^i{k,B)  =  e  jxf{y,k,B)ri-fxiiy,k,B)dfju}{x,y,k)ui{x,y,k)dn  (3) 

where  dQ  is  a  volume  element,  fj  is  the  unit  vector  along  the  direction  of  the 
incident  photon  polarization,  f  =  xSx  +  ydy  is  the  two-dimensional  radius 
vector,  and  subscripts  i,  f  stand  for  initial  and  final  states  respectively.  Now, 
if  we  assume  that  the  incident  light  is  polarized  along  the  y-direction  so  that 
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T]  =  dy,  the  above  equation  simplifies  to 

fW/2 

df,i{k,  B)  =  e<  Xf\y\Xi  >=e  yxf{y,  k,  B)xi{y,  k,  B)dy,  (4) 

J-WI2 

where  W  is  the  width  of  the  quantum  wire  along  the  y-direction.  One  should 
note  here,  that  if  there  is  no  magnetic  (or  electric)  field  applied,  the  envelope 
functions  Xi  ^ire  just  particle-in-box  states  and  the  dipole  moment  in  Eq.  (4) 
is  non  zero  only  for  the  transitions  between  subband  states  of  opposite  parity. 
However,  this  is  obviously  not  the  case  when  a  magnetic  field  is  present.  It 
is  clear  from  Eq.  (4)  that  to  calculate  the  dipole  moments  in  the  presence  of 
a  magnetic  field,  all  we  need  to  compute  are  the  wave  functions  Xf,i{y>  •®) 
at  a  given  magnetic  field  B,  for  given  magnetoelectric  subbands  /  and  i,  and 
for  a  given  wave  vector  k.  This  is  achieved  via  a  numerical  (finite  difference) 
solution  of  the  Schrodinger  equation  for  the  y-component  of  the  wavefunction 

-  (f )  M  +  ^^kx{y)  -  k'^xiy)  =  0  (5) 

with  I  being  the  magnetic  length  given  by  I  =  yJh/eB,  assuming  hardwall 
boundary  conditions 

xiy  =  W/2)  =  xiy  =  -W/2)  =  0  (6) 

and  following  the  prescription  of  Ref.  [8].  Once  this  is  done,  we  can  calculate 
the  dipole  moment  in  Eq.  (4)  for  any  chosen  intersubband  transition,  at  any 
chosen  magnetic  field  and  for  any  chosen  wave  vector. 

The  absorption  of  both  the  fundamental  firequency  (pump)  and  its  second- 
order  harmonic  is  very  important  when  considering  frequency  conversion  with 
low  insertion  loss.  In  general,  it  is  desirable  to  have  large  absorption  coef¬ 
ficient  q{u)  for  the  pump  frequency  and  small  01(20;)  for  the  converted  fre¬ 
quency  so  that  the  latter  is  not  re-absorbed  to  cause  large  insertion  loss.  In 
order  to  obtain  the  absorption  coefficients  in  the  whole  range  of  frequencies 
and  for  different  values  of  a  magnetic  flux  density,  we  need  to  calculate  the 
first-order  susceptibility  as  follows 


where  we  have  used  the  same  notation  as  in  Eq.  (2).  The  imaginary  part  of 
is  related  to  the  absorption  coefficient  while  the  real  part  is  related 
to  the  refractive  index. 

II.  Results 

We  now  present  results  of  our  calculations.  The  physical  parameters  used 
for  the  numerical  calculations  are  relevant  for  a  GaAs  quantum  wire  with 
relative  dielectric  constant  =  12.9,  and  effective  masses  =  0.067mo  and 
rrifi  —  0.5 Wo  where  ruo  is  the  free  electron  mass. 


Figure  1:  The  dipoles  of  three  inter-subband  transitions  as  functions  of  the 
applied  magnetic  field.  The  induced  dipole  dei-ez  peaks  at  a  magnetic  flux 
density  of  5.3  tesla. 

Fig.  1  presents  the  dipole  moments  for  the  lowest  intraband  transitions  as 
a  function  of  magnetic  flux  density.  At  zero  magnetic  field,  a  non-vanishing 
dipole  matrix  element  occurs  only  for  transitions  between  states  of  opposite 
parity  (el-e2,  e2-e3)  as  expected  from  Eq.  (4).  Transition  dipole  des-ei 
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has  a  non-monotonic  dependence  on  the  magnetic  field.  This  transition  is 
forbidden  at  zero  field  since  the  wave  functions  of  the  first  and  third  subband 
have  the  same  parity.  At  low  and  moderate  magnetic  fields,  the  parities  are 
altered  by  the  skewing  of  the  wavefunctions.  The  skewing  effect  of  the  wave 
functions  of  the  first  and  the  third  subbands  is  shown  in  Fig.  2  (top).  Its 
degree  depends  on  a  subband  number,  which  causes  a  breaking  of  inverion 

symmetry  and,  consequently,  non-zero  value  of  dipole  matrix  element  des-ei 
for  otherwise  forbidden  transition. 


Figure  2:  (Top  left  and  right).  Skewing  of  the  wave  functions  of  the  first 
and  third  subbands  in  a  magnetic  field.  The  left  panel  corresponds  to  zero 
magnetic  flux  density  and  the  wave  functions  are  particle-in-a-box  states. 
The  right  panel  corresponds  to  a  flux  density  of  3  tesla  and  the  wave  functions 
are  those  of  ’’edge  states”.  (Bottom  left  and  right).  The  energy  spacing 
between  the  mth  and  nth  subbands  vs.  wave  vector  A:  at  a  magnetic 
fl^  density  B=1  tesla  (left)  and  B=3  tesla  (right).  The  lowest  curve  (at 

k— 0)  corresponds  to  el-e2,  the  intermediate  curve  to  e2-e3,  and  the  highest 
to  el~e3.  ° 


The  dipole  moment  reaches  a  maximum  of  about  6  e-A  and  then  de- 
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creases.  This  later  decrease  is  related  to  the  following  effect.  For  some  fixed 
wave  vector  k  (electron  velocity),  a  sufficient  increase  in  the  flux  density  B 
forces  the  traversing  states  (“skipping  orbits”  or  “edge  states”)  to  condense 
into  closed  cyclotron  orbits  (Landau  levels)  which  are  no  longer  skewed  by 
the  magnetic  field  to  the  wire  edge  since  they  have  no  translational  velocity 
and  hence  experience  no  Lorentz  force.  While  edge  states  have  a  skewed 
wave  function  which  is  not  symmetric  in  space,  cyclotron  orbits  have  a  wave 
function  that  is  symmetric  about  the  orbit  center.  Note  that  the  orbit  cen¬ 
ter’s  coordinates  depend  only  on  k  and  B.  Therefore,  at  a  fixed  k,  the  wave 
functions  of  the  first  and  third  Landau  levels  are  symmetric  about  a  common 
center.  Whenever  this  kind  of  symmetry  holds,  dgs-ei  vanishes.  Therefore, 
the  dipole  moment  des-ei  decreases  gradually  to  zero  at  high  magnetic  field 
with  the  onset  of  Landau  condensation  (5  tesla  for  this  wire  dimensions) . 

In  Fig.  3,  we  plot  the  absolute  values  of  as  a  function  of  photon 
energy  for  two  different  values  of  the  magnetic  field.  In  our  numerical  cal¬ 
culations  we  have  used  a  dilute  carrier  concentration  of  ^=10^^  cm~^  which 
allows  us  to  neglect  high  density  effects  such  as  screening  and  bandgap  renor¬ 
malization.  Both  susceptibility  curves  have  pronounced  three-peak  resonant 
structure  which  corresponds  to  two  one-photon  transitions  el-e2  (at  72  meV) 
and  e2-e3  (at  124  meV)  and  one  two-photon  transition  el-e3  (at  100  meV) 
between  magneto-electric  subbands.  These  three  peaks  have  different  broad- 
enings  because  the  sum  in  Eq.  (2)  represents  an  integral  effect  of  all  direct 
transitions  with  different  values  of  electron  wave  vector  k  and  because  of  the 
complex  dependence  of  the  subband  spacing  on  k  and  B  (see  Fig.  2 
(bottom)).  The  latter  also  gives  rise  to  an  uneveness  in  the  second  order 
susceptibility  peaks.  The  peak  value  of  the  second  order  susceptibility  is 
^(2)  =  14.5A/V’  for  1  tesla  field  (left  panel);  and  =  43.1A/V’  for  3  tesla 
(right  panel).  For  comparison,  the  nonlinear  susceptibility  of  electric  field 
biased  GaAs  quantum  wells  (W=92  A)  -  calculated  theoretically  and  mea¬ 
sured  experimentally  in  Ref.  [9]  -  was  x^^^=240A/V  for  an  electric  field  of  36 
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Figure  3:  Second  order  susceptibility  as  a  function  of  the  photon  energy 
for  two  values  of  the  biasing  magnetic  field.  The  left  panel  corresponds  to  a 
magnetic  flux  density  B=1  tesla,  the  right  panel  to  B=3  tesla.  The  maximum 
values  of  the  curves  are  li.bA/V  for  1  tesla  field  and  AZX/V  for  3  tesla 
field. 

kV/cm.  The  carrier  concentration  used  in  their  calculations  was  N=  5x10^^ 
cm“^.  Adjusted  to  that  carrier  concentration,  the  second  order  susceptibility 
for  a  3-tesla  magnetic  field  is  about  215A/V’  compared  to  240A/V’  of  Ref. 
[9].  This  shows  that  relatively  weak  magnetic  fields  in  quantum  wires  can 
produce  similar  magnitudes  of  as  rather  strong  electric  fields  in  quantum 
wells. 

Fig.  4  shows  the  dependence  of  the  imaginary  part  of  the  first-order  sus¬ 
ceptibility  Im{x^^^)  as  a  function  of  photon  energy  for  two  different  values  of 
the  magnetic  field.  We  have  used  the  same  wire  dimensions  and  carrier  con¬ 
centrations  for  this  plot  as  in  the  previous  one.  The  same  physics  pertinent 
to  the  previous  plot  explains  different  broadening;  however,  the  peaks  are 
now  attenuated  because  of  averaging  over  different  transition  probabilities. 
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Figure  4;  Imaginary  part  of  the  first-order  susceptibility  as  a  function  of  the 
photon  energy  for  B=1  tesla  (left  panel)  and  B=3  tesla  (right  panel). 

Since  is  related  to  the  absorption  coefficient  a(u!)  as 

a{(j)  =  (8) 

one  can  estimate  the  absorption  over  the  whole  frequency  range.  At  res¬ 
onant  photon  energies  of  72  meV  and  124  meV,  the  absorption  coefficient 
a  =  1.5  10^  cm~^  and  a  =  4.5  10^  cm~^,  respectively.  It  is  clear  from  the 
figure,  that  the  absorption  coefficient  at  twice  these  frequencies,  Q!(2o;)  is 
much  less.  This  implies  that  a  large  portion  of  the  pump  energy  at  these 
resonant  frequencies  will  be  absorbed  by  the  structure  and  converted  into 
second  harmonic  signal  which  will  not  be  significantly  re-absorbed. 

Another  important  factor  for  efficient  second  harmonic  generation  is  phase 
matching.  Since  the  refractive  index  n(u>)  of  most  materials  is  frequency  de¬ 
pendent,  the  following  inequality  holds  n(a;)  n(2a;).  As  a  result,  the  coher¬ 

ence  length  Ieoh  =  \j/4(n2u  —  n(j)  for  GaAs  (typical  non-birefringent  crystal) 
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varies  between  10;im  and  lOO^im.  The  efficiency  of  nonphase-matched  inter¬ 
actions  are  about  10“®  times  less  than  that  of  the  phase-matched  interactions 
over  a  length  scale  of  1  cm. 


Figure  5:  Real  part  of  the  first-order  susceptibility  as  a  function  of  the  photon 
energy  for  B=1  tesla  (left  panel)  and  B=3  tesla  (right  panel). 

Using  a  magnetic  field  as  an  additional  degree  of  freedom,  we  may  try 
to  adjust  n{u).  In  Fig.  5  we  present  the  dependence  of  the  real  part  of 
the  first-order  susceptibility  i?e(x^^^)  as  a  function  of  photon  energy  for  two 
different  values  of  the  magnetic  field.  Using  the  relation 

An{u)  =  27rRe{x^^^),  (9) 

and  Eq.  (8),  one  can  determine  the  frequencies  where  two  conditions  simul¬ 
taneously  hold:  n{u)  «  n{2u)  and  Q:(a;)  »  Q!(2a;).  For  the  3-tesla  field,  this 
frequency  hoj  corresponds  to  75  meV. 
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IV.  Conclusion 


We  have  theoretically  studied  second  harmonic  generation  in  a  semicon¬ 
ductor  quantum  wire  biased  with  a  magnetic  field.  A  strong  second-harmonic 
component  of  the  dielectric  susceptibility,  due  to  the  dipoles  associated  with 
transitions  between  magneto-electric  subbands,  is  found.  We  have  also  calcu¬ 
lated  absorption  coefficient  and  refractive  index  in  the  appropriate  frequency 
range  to  assess  the  efficiency  of  frequency  conversion  and  insertion  loss.  We 
have  shown  that  a  magnetic  field  can  be  used  as  an  additional  degree  of 
freedom  in  optimizing  second  harmonic  generation  efficiency.  This  may  have 
important  applications  in  nonlinear  optics. 
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